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OPERATING  CHARACTERISTICS  OF  CROSSCORRELATOR 
WITH  OR  WITHOUT  SAMPLE  MEAN  REMOVAL 

INTRODUCTION 

The  detection  of  weak  signals  in  two  channels  is  often  accomplished  by 
crosscorrel ating  the  two  waveforms  and  comparing  with  a  threshold.  For  the 
case  where  a  large  number  of  independent  products  are  added  to  yield  the 
correlator  output,  the  central  limit  theorem  is  often  employed,  with 
questionable  validity  for  low  false  alarm  probabilities,  i.e.  large  thresholds. 
Also,  this  approximation  may  not  be  valid  for  intermediate  numbers  of  terms 
added. 

Here  we  wish  to  get  exact  operating  characteristics  for  the 
crosscorrelator,  namely  detection  probability  vs.  false  alarm  probability, 
even  for  probabilities  as  low  as  IE— 10 .  In  particular,  we  desire  results  for 
an  arbitrary  number  of  products  summed,  for  any  degree  of  correlation  between 
corresponding  individual  samples  of  the  two  channel  inputs,  and  for  any  input 
signal  and  noise  power  levels. 

Furthermore,  it  sometimes  happens  that  the  two  input  channels  contain  dc 
components,  which  can  be  considered  either  desirable  or  otherwise,  depending 
on  the  application.  Here  we  will  consider  these  dc  components  as  nuisance 
terms  and  will  subtract  them  out  prior  to  crosscorrelation.  More  precisely, 
since  the  actual  values  of  the  dc  components  in  each  channel  will  generally  be 
unknown,  we  will  estimate  them  via  the  sample  means  (over  the  available  record 
lengths)  and  subtract  these  estimates  from  the  available  data.  This 
subtraction  feature  creates  new  random  variables,  all  of  which  are 
statistically  dependent  on  each  other,  and  thereby  significantly  complicates 
the  analysis.  Nevertheless,  this  crosscorrelation  of  the  sample  ac  components 
of  the  input  channels  is  encountered  in  practical  situations,  and  in  one  recent 
study  [1],  it  was  in  fact  the  generalized  likelihood  ratio  detector  under  two 
different  realistic  scenarios.  Accordingly,  it  merits  study  and  accurate 
quantitative  evaluation  of  performance  capability. 
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More  generally,  we  consider  subtraction  of  scaled  versions  of  the  sample 
means  of  each  channel  prior  to  multiplication  and  summation.  Then  as  special 
cases,  we  can  investigate  the  crosscorrelator  with  or  without  sample  mean 
removal,  or  any  intermediate  case  of  interest. 

The  major  analytical  result  here  is  a  closed  form  for  the  characteristic 
function  of  the  correlator  output,  in  the  most  compact  form  involving  only  two 
rooting  operations  and  one  exponential.  Although  this  processor  could  be 
analyzed  by  the  general  method  given  in  [2],  in  terms  of  the  eigenvalues  and 
eigenvectors  of  a  correlation  matrix,  it  would  be  less  accurate  and 
considerably  more  time  consuming,  even  with  computer  aid,  especially  for  a 
large  number  of  terms  summed.  The  actual  numerical  procedure  adopted  here  for 
proceeding  from  the  characteristic  function  to  the  exceedance  distribution 
functions  (false  alarm  and  detection  probabilities)  is  that  given  in  [3],  and 
utilized  to  advantage  in  [2,3,4]. 
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PROBLEM  DEFINITION 


INPUT  STATISTICS 


The  two  channel  inputs  to  the  crosscorrelator  are  synchronously  sampled 
in  time,  yielding  random  variables  fun]  ^  and  (v^,  where  N  is 
the  total  number  of  data  samples  taken  in  each  channel.  These  random 
variables  are  Gaussian  with  the  following  statistics: 


means 


un  "  V  vn  =  V 


_ p  2  _  2 

variances  (u  -u„)  =  a,  (v  -v„) 

'  n  n'  u  v  n  n' 


2 

V 


covariances  (un-un) (vn-vn)  =  p auav. 


all 

independent  (1) 

of  n. 


(An  overbar  denotes  a  statistical  average.)  That  is,  the  means  and  variances 
in  each  channel,  although  different,  do  not  change  with  time,  and  the  degree 
of  correlation  between  channels  is  constant.  Also 


um  is  statistically  independent  of  un  if  m^n, 

vm  is  statistically  independent  of  vn  if  m^n,  (2) 

um  is  statistically  independent  of  vn  if  m^n. 

However,  un  and  vn  are  statistically  dependent  on  each  other,  for  all  n, 

to  the  extent  p  indicated  in  (1). 

A  SIGNAL  AND  NOISE  MODEL 


To  better  fix  the  mathematical  definitions  above,  consider  in  this 
subsection  the  following  possible  signal  and  noise  model: 


un  =  uu  +  +  ud(n) 


vn  =  uv  +  vs(n)  +  vd(n) 


for  1  <  n  <  N, 


(3) 
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where  random  signal  components  u$(n),  v$(n)  are  zero-mean  and  partially 
correlated  with  each  other: 


for  all  n.  (4) 


Thus  Su,  sy  are  the  powers  of  the  random  signal  components  in  each  channel. 

Also,  the  random  noise  disturbances  ud(n),  vd(n)  in  (3)  are  zero-mean 
and  independent  of  each  other: 


for  all  n.  (5) 


Thus  Du,  Dy  are  the  powers  of  the  random  noise  disturbances  in  each 
channel.  Finally,  except  for  the  statistical  dependencies  indicated  in  (4) 
between  us(n)  and  v$(n),  all  the  4N  random  components  in  (3)  are 
independent  of  each  other. 

For  this  particular  signal  and  noise  model  in  (3)— (5) ,  the  master 
parameters  in  (1)  take  the  special  form 


(6) 


from  which  there  follows 


(7) 
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where  the  signal-to-noise  ratios  (per  sample)  of  the  random  components  in  (3) 
have  been  defined  as 


R  =  TT 
u  D. 


^(n) 

ud(n) 


for  all  n. 


(8) 


Thus  the  parameters  ou,  0  p  in  (1)  depend  only  on  the  statistics  of  the 
random  components  in  model  (3),  and  not  on  the  dc  components  pu  and  yv. 
Observe  that  even  if  ps=1  and  Ru=oo,  p  would  still  be  less  than  1;  the  one 
noisy  channel  prevents  full  correlation  between  inputs. 


CROSSCORRELATOR  OUTPUT 


We  define  the  sample  ac  components  of  each  channel  of  the  crosscorrel ator 
as 

N 

“n  '  un  -  *  JE  \ 
m=l 

N 

’n  -  vn  -  S  2  vm 
m=l 

where  we  have  subtracted  the  corresponding  sample  means  from  each  and  every 
data  sample.  Thus  ^ul ^  and  {v^}^  have  zero-means  and  have  statistics 
completely  independent  of  the  unknown  actual  values  of  input  means  pu,  p^. 
However,  in  trade,  we  now  must  deal  with  a  new  set  of  2N  random  variables,  all 
of  which  are  statistically  dependent  on  each  other;  this  is  the  feature  which 
complicates  the  ensuing  analysis.  The  test  statistic  (decision  variable)  of 
interest  is  the  crosscorrelator  output  after  sample  mean  removal, 

N  N  N  N 

"  "2_Vn  -  2  Vn  uml>  vn  •  (10) 

n=l  n=l  m=l  n=l 
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which  is  independent  of  the  actual  unknown  values  of  input  means  uu  and  nv> 
If  we  knew  the  input  means,  we  could  subtract  them  directly  and  not  have  to 
resort  to  sample  means. 


More  generally,  we  consider  the  modified  channel  components 


N 


a  % 

U„  =  U„  -  TT  U 

n  n  N  n 


m=l 

N 


for  1  <  n  <  N 


;n  '  vn  -  t  vn 
m=l 


(ID 


and  the  crosscorrelator  output 


N  N 


q  -  ]>  Vn  -  2  Vn  -  N  ^  ura  ^  vn  • 

n=l  n=l  m=l  n=l 


(12) 


instead  of  (9)  and  (10).  Scale  factors  o.  and/or  s  in  (11)  may  be  unequal  to  1; 
the  final  parameter  y  in  (12)  is  given  by 


y  =  a  +  8-ct&  =  l-  (a-l)(p-l). 


(13) 


The  case  of  y=0  in  (12)  obviously  corresponds  to  the  case  of  no  sample  mean 
removal.  On  the  other  hand,  if  either*  a=l  or  b=1,  then  y=l»  and  we  have 
removal  of  the  sample  mean;  i.e.,  (12)  reduces  to  (10).  We  shall  be  interested 
here  in  the  analysis  of  the  general  case  represented  by  (12),  for  arbitrary  y 


*  It  is  demonstrated  in  appendix  A  that  if  scale  factor  a=l  but  b^I, 
correlator  output  q  is  completely  independent  of  v  g. 
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CHARACTERISTIC  FUNCTION  OF  CROSSCORRELATOR  OUTPUT 


DERIVATION 

We  express  the  collection  of  random  variables  in  (1)  and  (2)  in  column 
matrix  form  according  to 


T  T 

U  =  [u1  u2  ...  uN]  ,  V  =  [v1  v2  ...  vN]  ,  (14) 

where  superscript  T  denotes  transpose.  The  crosscorrelator  output  q  in  (12) 
can  then  be  written  as  quadratic  form 

q  =  UTQ  V  ,  (15) 

where:  NxN  matrix 

Q  =  I  -  J  1  1T  ,  (16) 

I  is  the  NxN  identity  matrix,  and 

1  =  [1  1  ...  if  (17) 

is  a  Nxl  column  matrix  of  ones. 


Since  U  and  V  are  Gaussian,  their  joint  probability  density  function  is, 
in  terms  of  the  parameters  in  (1), 


P(U,V) 


2(1-p2) 


★ 


*^(U  -  uul)T(U 
°u 


uui)  +  Vv  -  \1)T(V  -  V0 

<7 


^(u  -  »U1,T(V 


“v1 


The  characteristic  function  of  correlator  output  q  in  (15)  is  then  given  by 
the  statistical  average 
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fq(f )  =  exp(i?q)  =  exp(i?UTQV)  = 


=  jJdU  dV  p(U,V)  exp(if llTQV)  = 


-r]> 


dV  exp 


if  UTQV - * 


2(1-7) 


*<  Vu  -  »u1)T(u  -  UU1}  +  Vv  -  N1)T(v  -  N1}  -  ^T(u  -  \,1)T(V  -  Ni: 

au  av  U  V 


.(19) 


At  this  point,  in  order  to  evaluate  this  2N-fold  integral,  we  employ  the 
general  integral  result  (B-2)  and  (B-6)  in  appendix  B,  identifying  the 
matrices  there  as 


A  =  I,  B  =  0  1  0  I,  C  =  if Q  + 

W1  <-\  <-  /  i 

au(l-P  ) 


%( IV) 


a  a  (1-p^) 
U  Vv  ' 


I, 


a  u  —  pc  u  a  u  -  pc  u 

D  _  VMU  M  UMV  ^  E  u  V  V  U  ^ 


2  / 1  2  \ 
auav(!-p  ) 


2n  2% 
a  c  ( i-p  ) 
u  vv  ' 


(20) 


We  also  need  the  following  auxiliary  results  for  special  matrix  forms;  namely, 
for  arbitrary  scalars  c^,  C2,  the  matrix  determinant 

det(c1I  +  c2  1  1T)  =  ci_1(ci  +  Nc2 ) ,  (21) 

and  the  matrix  inverse 


(cxI 


_1_  i  T  C2  i  XT 

cl  c1(c1  +  Nc2) 


(22) 


Employment  of  appendix  B  and  ( 20 ) — ( 22 )  then  yields,  after  a  very  considerable 
amount  of  effort,  a  closed  form  for  the  characteristic  function  in  (19)  (in 
its  most  compact  form) 
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exp  if 


Gj_  +  ifG 


1  p  2 


(23) 


N=T 


1 

2 


}  -  *Ei 


where 


Fj_  =  Ejd-y),  F2  =  E2(1-t)?  , 


(24) 


The  square  roots  in  (23)  are  principal  value,  being  +1  at  £=0.  This 
characteristic  function  has  four  branch  points  and  two  essential  singularities 
which  overlap  two  of  the  branch  points;  the  complexity  of  this  characteristic 
function  of  q  precludes  tractable  analytical  results  for  the  probability 
density  function  or  exceedance  distribution  function  of  the  correlator  output, 
except  in  very  special  cases.  Nevertheless,  since  the  characteristic  function 
in  (23)  is  easily  numerically  evaluated  with  computer  aid,  it  readily  lends 
itself  to  the  procedure  presented  in  [2,3].  A  program  for  the  evaluation  of 
the  cumulative  and  exceedance  distribution  functions  corresponding  to 
characteristic  function  (23)-(24)  is  given  in  appendix  C  for  arbitrary  values 
of 


N,  number  of  terms  summed 
y ,  scale  factor  in  sample  mean  removal 

uu,  mean  in  u-channel 

uv,  mean  in  v-channel 

<*u,  standard  deviation  in  u-channel 
o  ,  standard  deviation  in  v-channel 
p,  correlation  coefficient  between  channels. 
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A  sample  plot  of  the  cumulative  and  exceedance  distribution  functions  for  a 
typical  selection  of  numerical  values  for  the  above  parameters  is  also 
presented  in  appendix  C. 

CUMULANTS  OF  CORRELATOR  OUTPUT 

By  taking  the  natural  logarithm  of  the  characteristic  function  in  (23) 
and  expanding  in  a  power  series  in  J ,  the  cumulants  of  random  variable  q  can 
be  extracted: 


(25) 


where  here 


(26) 


S=p+1,  D  =  p  -  1  . 


In  particular,  the  mean  and  variance  of  q  are  available  by  using  n=l  and  2 
respectively  in  (25): 


uq  =  (N-Y)pau°v  +  n(1“y)uuuv> 


(27) 


SPECIAL  CASE  OF  Y=l,  SAMPLE  MEAN  REMOVAL 

For  y=1>  the  general  characteristic  function  in  (23)  reduces  to 


N-l 
-  T 
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where  E^  and  E2  are  still  given  by  (24),  and  are  independent  of  means  uu 
and  u  as  shown  earlier.  The  cumulants  in  (25)  reduce  to 

*q(")  =^(n-l)l(N-l)(auav)n[(p+l)n  +  (p-l)n]  ,  (29) 

and  in  particular,  the  mean  and  variance  of  q  are 

yq  =  (N-l)pauV 

°q  =  (N-l) (1+p2)°uctv  •  (30) 

SPECIAL  CASE  OF  Y=0,  SAMPLE  MEAN  NOT  REMOVED 


For  Y=0,  the  characteristic  function  in  (23)  reduces  to 


/  2  \~N/2 
fq (F;y=0)  =  (1  -  IJEj  +  *%)  exp 


g{0)  +  iSG^o) 

15- - f- 


(31) 


where 


E1  “  2pauV  E2  =  Vv(1"p2)’ 


.(0)  ..  r(o)  L.,  2  2  2  2  0  . 

’1  =  NVv’  G2  =  ?"(Vv  Vu  "  2pauCTvuuuv) 


(32) 


The  cumulants  are  obtained  by  setting  Y=0  in  (25),  and  in  particular,  the  mean 
and  variance  of  correlator  output  q  are 


=  N(poL1ov  +  uuiiv). 


-i 


N l(l+p  ) 


2v  2  2 


a  a 
U  V 


2  2 
(J  u 
ITV 


2  2 
a  u 
vMu 


2pau°vpupv]< 


(33) 
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INTERRELATIONSHIP  OF  TWO  SPECIAL  CASES 


Let  the  general  characteristic  function  in  (23)  be  denoted  by 
fq(|;  N,  y,  yu,  yv).  We  have  already  seen  the  expression  for 
fq(f;  N,  1,  uu,  pv)  in  (28).  At  the  same  time,  from  (23)  and  (24),  there 
follows 


fq(|;  n-1,  o,  o,  o. 


(34) 


which  is  identical  to  (28).  That  is. 


fq(J;  N,  1,  v  pv)  =  fq(?;  N-1,  0,  0,  0). 


(35) 


Thus  the  character! Stic  functions  of  the  two  following  random  variables  are 
identical : 


(1)  Sum  of  N  terms  with  sample  mean  removal,  and  the  true  means 
arbitrary, 

(2)  Sum  of  N-1  terms  without  sample  mean  removal,  but  the  true  means 

zero.  (36) 


The  removal  of  the  sample  means  has  eliminated  the  dependence  of  the  correlator 
output  on  the  unknown  means  but  has  reduced  the  number  of  degrees  of  freedom  by  1 

SPECIALIZATION  TO  THE  SIGNAL  AND  NOISE  MODEL 

For  general  scaling  factor  y  and  arbitrary  input  means  pu,  pv,  and  for 
the  model  introduced  earlier  in  (3)-(6),  the  general  characteristic  function  of 
the  correlator  output  is  still  given  by  (23),  but  with  the  parameters  in  (24) 
now  specialized  to  the  form 
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E,  =  2p  (S  S  )1/2,  E„  =  D  D  +D  S  +D  S  +S  S  (1-p2), 

1  s'  u  v'  *2  u  v  u  v  v  u  u  v'  ps'’ 

Fx  =  E^I-y),  F2  =  E2(1-y)2, 

Gj  -  N(l~,)Vv,  G2  .  |n(1-v)2[(Su*Ou)4*(S/Dv)4-2»s(SuSv)1,2uuI.v]  .  (37) 

The  general  n-th  cumulant  is  still  given  by  (25);  however,  the  use  of  (6) 
allows  for  determination  in  terms  of  the  fundamental  quantities  of  the  signal 
and  noise  model,  namely  Su,  Sy,  Du,  Dy,  P$  defined  in  (4)-(5).  In 
particular,  the  mean  and  variance  of  correlator  output  q  are 

uq  -  (N-y)»s(SuSv)1/2  +  N(1-y)uuuv, 

"q  -  (N-2rY2)[0uDv*0uSvt0vSu*(l*»2)SuSv]+ 

*  M(1-y)2[(S/Du)u^(S/Dv)„^2»s(SuSv)1/2Vv].  (38) 
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ANALYTIC  RESULTS  FOR  Y=l>  SAMPLE  MEAN  REMOVAL 

In  this  section  and  the  next,  we  will  confine  attention  solely  to  the 
case  of  scale  factor  y=l.  The  characteristic  function  of  the  crosscorrelator 
output  q  follows  from  (28)  and  (24)  as 

N-l  N-l 

2  “  2 
fq(?)  =  [1  -  1JEj_  +  ?%]  =  [1  -  152Pauav  +  ?2a2a2(l-p2)] 

_  M. 

=  {[l+'iFouav(l-p )]  [l-i?®uov(l+p)j  2  for  y=l  >  (39) 

where  we  must  have  N^2.  We  observe,  for  later  numerical  use  in  appendix  D, 
that  since  |l  ±  ifbj  =  (1  +  i$  monotonical ly  increasing  for  J  >0, 

then  Jf  (J)|  is  monotonical ly  decreasing  for  all  J>0  and  any  N,  cu,  oy,  p. 

GENERAL  PROBABILITY  RESULTS 

The  cumulants  of  q  have  already  been  listed  in  (29)  and  (30).  The 
probability  density  function  corresponding  to  characteristic  function  (39)  is 
given  by  [5,  6.699  12] 


where  K  (z)  is  a  modified  Bessel  function  of  the  second  kind  [6,  section  9.6]. 
If  the  number  of  terms  added,  N,  to  yield  correlator  output  q,  is  odd,  simple 
relations  for  the  probability  density  function  in  (40)  can  be  obtained 

[6,  10.2.15  and  10.1.9,  last  equation];  letting  n  =  for  N  odd,  we  find  the 

exact  result 


15 


TR  7045 


,  x  (l-p^)n  /  pu  —  lu | 

pq(u)  =  r~  ,n  -  exp  (  — 7T2 


2au%  4"  nl  V<Vv(1"p  ]/  m=0 


( 2n— m ) 

(n-m)!  ml 


2  lu] 


vm 


auav(l-P  V 


for  all  u;  n  =  — ,  N  =  3,  5,  7,  ...  . 

For  example,  for  N=3,  we  have  n=0,  yielding 

p  (u)  =  ^ —  exp  \  for  an  u< 


(41) 


(42) 


The  corresponding  cumulative  distribution  function  for  N=3  is 


Pq(u) 


-  f 

—  oo 


dt  pq(t) 


1-p 


exp  | 


,Vv(1-p), 


for  u<0. 


(43A) 


while  the  exceedance  distribution  function  is 


1  - 


Pq( 


+  00 

u)  .[d 


dt  pq(t) 


1+p 


exp 


-  u 


,auav(1+p) 


for  u>0. 


(43B) 


This  dichotomy,  of  presenting  the  cumulative  distribution  function  for 
negative  arguments,  and  the  exceedance  distribution  function  for  positive 
arguments,  turns  out  to  be  notationally  convenient  and  physically  meaningful 
and  will  be  adopted  throughout  this  report. 


Although  closed  form  expressions  for  the  exceedance  distribution  function 
corresponding  to  probability  density  function  (40)  are  not  available  for 
general  N,  the  use  of  [6,  9.7.2]  on  (40)  leads  to  the  dominant  term  in  the 
asymptotic  expansion  of  the  exceedance  distribution  function: 


as  u  +o6  . 


For  N=3,  this  is  precise;  see  (43B). 


(44) 
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POSSIBLE  NORMALIZATIONS  OF  q 

If  we  define  a  normalized  random  variable 


x  =  S 


pl/2  ‘  o  1/2  * 

2  Vv(1-p  ) 


(45) 


then  the  characteristic  function  of  x  is  given  by  (39)  as 


1/2, 


fx(F)  =  yi/Ep)  = 


1-ifEj/E^2  +  f2] 


N-l 


52p(lV) 


-1/2 


,<] 


N-l 

~r 


which  has  only  two  fundamental  parameters,  namely,  N  and  p, 


(46) 


A  second  possibility  is  the  random  variable  defined  by 


y  = 


a  a 
U  V 


for  which  characteristic  function 


1  -  i?2p 


N-l 

2 


(47) 


(48) 


also  depends  only  on  N  and  p.  However,  neither  of  the  normalizations,  (45)  and 
(47),  are  of  interest  to  us  here;  an  alternative  normalization  and  reasons  for 
its  selection  are  given  below. 
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SPECIALIZATION  TO  THE  SIGNAL  AND  NOISE  MODEL 

For  the  model  presented  earlier  in  (3)-(8),  the  original  E^, 
parameters  in  (24)  take  the  form  already  given  in  the  upper  line  of  (37). 
Let  a  normalized  random  variable,  relative  to  the  additive  random  noise 
disturbances,  be  defined  according  to 


h  = 


<W 


1/2  ’ 


(49) 


see  (5).  This  normalization  for  the  particular  signal  model  (3)  is  different 
from  both  x  and  y  in  the  general  case  above.  The  reason  we  employ  h  is  that 
the  normalization  depends  only  on  the  power  of  the  additive  noise 
disturbances,  and  not  on  the  signal  strengths  or  correlation  coefficients; 
this  is  consistent  with  a  system  which  monitors  the  noise-only  background  and 
sets  a  threshold  for  a  desired  false  alarm  probability. 


The  characteristic  function  of  the  normalized  random  variable  h  in  (49) 
is  given  by 


N-l 


V?>  -  fq(?'<W1/2)  ‘ 


1  -  if 


<W 


T7?+£ 


2  l2 

D  D 
u  v 


_  N-l 

=  [l  -  i J2a  +52(0?-a2^]  2 

where  we  define  auxiliary  parameters  here  as 


°  “  ps(RuRv)1/2  ’  8  =  C(l+Ru)(l+Rv)]1/2  • 


(50) 


(51) 


Here  we  used  (39),  (37),  and  (8).  This  characteristic  function  in  (50) 
depends  on  the  four  fundamental  parameters  N,  p$,  Ru,  R  where  the 
latter  two  quantities  are  the  signal-to-noise  ratios  per  sample  of  the  random 
components  of  model  (3);  see  (8). 
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Reference  to  (40)  reveals  that  the  probability  density  function  of  h 
corresponding  to  characteristic  function  (50)  is  given  by 


Ph(u) 


(s2-«¥'2] 


for  all  u. 


* 


(52) 


F.or  N  odd,  alternative  forms  are  available  from  (41),  if  desired.  The 
asymptotic  behavior  of  the  exceedance  distribution  function  of  h  follows  in  a 
manner  similar  to  that  used  for  (44): 


N-3 


The  cumulants  of  h  follow  from  (29),  (26),  and  (6)-(8): 

Xh(n)  =  -^(n-l )  1  (N-l )[ (o+8 )n  +  (a-8)n]  ,  (54) 

and  in  particular,  the  mean  and  variance  of  h  are 

uh  =  (N-l)a  =  (N-1)ps(RuRv)1/2  , 

a2  =  (N-l)(a2+82)  =  (N-l)[l  +  Ru+  Rv+  R^U+p2)]  .  (55) 

The  two  parameters,  a  and  e,  are  given  here  by  (51),  in  terms  of  the 
fundamental  quantities  Ru,  R  p$  of  the  signal  and  noise  model. 

REDUCTION  TO  IDENTICAL  SIGNAL  COMPONENTS 

At  this  point,  we  will  further  specialize  the  results  for  the  signal  and 
noise  model  in  the  above  subsection.  We  presume  that 
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Ru  =  Rv  =  R  and  p$  =  1  , 


(56A) 


giving,  from  (51), 


a  =  R,  6  =  1+R;  (56B) 

that  is,  the  signal-to-noise  ratios  in  the  two  channels  are  equal,  and  the  two 
channel  signals  are  fully  correlated.  This  corresponds  physically  to  a  case 
where  the  random  signal  components  in  (3)  are  identical,  us(n)  =  vs(n), 
and  the  independent  random  noise  disturbances  have  the  same  power  level.  This 
situation  will  hold  for  the  rest  of  this  section  and  all  of  the  next  section 
where  the  graphical  results  are  presented. 

Equations  (50)  and  (56B)  then  yield  the  characteristic  function  for 
normalized  random  variable  h  in  (49)  as 


N-l 

2  ’  2 

fh(?)  =  [1  -  if 2R  +f£(  1+2R)] 

_  M 

=  [(1  +  1f)(l  -  ifd+2R))]  2  .  (57) 

The  cumulants  in  (54)  reduce  to 

^(n)  =  ^(n-l)I(N-l)[(l+2R)n+(-l)n]  ,  (58) 

and  in  particular,  the  mean  and  variance  become 

Mh  =  (N-l)R,  a2  =  (N-l) (1+2R+2R2)  .  (59) 

All  the  above  statistical  descriptions  depend  only  on  the  two  parameters  R, 
the  per-sample  signal-to-noise  ratio,  and  N,  the  number  of  terms  added  to 
yield  the  correlator  output. 
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The  probability  density  function  for  h  follows  from  (52)  as 


and  the  asymptotic  exceedance  distribution  function  from  (53): 


(60) 


N-3 


An  important  word  of  caution  must  be  mentioned  at  this  point:  when  N  is  large, 
(61)  is  inadequate  for  evaluating  small  false  alarm  and  detection 
probabilities,  since  the  succeeding  terms  in  the  asymptotic  expansion 
contribute  significantly.  For  example,  when  R=0,  the  maximum  value  of  the 
dominant  term  (61)  occurs  when  u  =  (N-3)/2  which,  for  N=128,  yields  false 
alarm  probability  3.86E-21,  a  value  far  below  those  of  interest.  Thus  (61) 
has  limited  applicability,  being  best  for  small  N;  in  fact,  the  first 
correction  term  to  (61)  yields  the  multiplicative  factor 

i  ^  1+2R  (N-3) (N+3+4R) 

1  m - an - •  (62) 

p 

It  indicates  that,  for  large  N,  u  must  be  of  the  order  of  N  in  order  for 
the  dominant  term  (61)  to  be  fairly  accurate. 

Although  for  N  odd,  an  alternative  closed  form  to  the  probability  density 
function  (60)  of  h  is  available  from  (41),  the  exceedance  distribution  function 
will  generate  a  double  sum  and  be  rather  cumbersome  for  large  N.  On  the  other 
hand,  the  characteristic  function  in  (57)  decays  rapidly  with  f  when  N  is  large 
and  yields  very  nicely  to  the  numerical  approach  given  in  [2,3].  The  only  diffi¬ 
cult  cases  are  in  fact  those  for  small  N;  accordingly,  some  analytic  results  for 
N  =  2,  3,  4,  5  will  now  be  presented,  based  on  characteristic  function  (57). 
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GENERAL  DISTRIBUTION  INTEGRALS 


Suppose  a  random  variable  y  has  characteristic  function  fy(£).  The 
cumulative  distribution  function  of  y  can  be  written  as  a  contour  integral 
[3,  (5)«  (6)] 


P  (u)  =  - 

y 


l 

i2* 


f 


exp(-iuf) 


for  all  u  , 


(63) 


where  C+  is  a  contour  along  the  real  axis  of  the  complex  J-plane,  with  an 
upward  indentation  at  the  origin  f=0,  to  avoid  the  pole  of  the  integrand  there. 


Similarly,  the  exceedance  distribution  function  of  random  variable  y  can 
be  expressed  as 


1  "  Vu)  =  J  df  —f~  exp(-iuf)  for  all  u  ,  (64) 

C_ 

where  C_  is  a  contour  along  the  real  f  axis,  with  a  downward  indentation  at 

?=o. 


For  u<0,  both  contours  can  be  closed  in  the  upper-half  ?-pl ane,  since  the 
exp  term  furnishes  rapid  decay  there.  Similarly,  for  u>0,  both  contours  can 
be  closed  in  the  lower-half  f-plane,  to  realize  exponential  decay  on  the 
circular  arcs  tending  to  infinity. 

DISTRIBUTIONS  FOR  N=2 

From  (57),  the  characteri Stic  function  of  normalized  correlator  output  h 
is 


fh(f)  =  C(l+i?)(l-i5(l+2R))]"1/2  ,  (65) 

and  the  probability  density  function  follows  from  (60)  as 
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*Vu)  ~  u(l+2R)1/2  SXPfe)  f°r  311  U  * 

There  is  no  closed  form  for  the  indefinite  integral  of  a  KQ  function;  see 
[6,  11.1.8  and  11.1.9].  Instead,  we  use  (65)  in  (63)  and  move  the  contour 
upwards  until  it  wraps  around  the  branch  point  at  £=i  and  extends  vertically 
from  there;  this  is  in  fact  the  steepest  descent  direction  for  the 
exponential.  The  contributions  of  the  small  and  large  circular  arcs  tend  to 
zero  as  the  radii  tend  to  zero  and  infinity,  respectively.  Under  a  change  of 
variable,  there  follows  the  cumulative  distribution  function  in  the  form 


P  fui  =  1  C  dt  exp[u(l+t2)] 
h  11  l  (l+t2)[l+(l+2R)(l+t2)]1/2 


for  u  <  0  . 


(67) 


This  is  a  useful  exact  result  for  several  reasons:  the  integrand  decays 
rapidly,  has  no  cusps,  and  involves  only  elementary  functions  which  are  easily 
computed;  also  the  integral  is  a  sum  of  positive  quantities  and  retains 
significance  even  for  large  |u). 


In  a  similar  fashion,  if  characteristic  function  (65)  is  substituted  in 
(64)  and  the  contour  moved  down  and  wrapped  around  the  branch  point  at 
|  =  -i/(l+2R)  and  along  the  vertical  steepest  descent  direction  for  the 
exponential,  the  exceedance  distribution  function  becomes,  upon  a  change  of 
variable. 


Ph^u ^  “  |j(1+2R)  exp Ci+2r) 


+  06 
r 


dt  exp(-uL  ) 


[l+(l+2R)t2][l+(l+2R)(l+t2)] 


Yj2  for  u  >  0.  (68) 


This  is  useful  for  the  same  reasons  given  above. 

There  is  one  closed  form  result  possible;  namely,  for  u=0,  direct 
integration  of  probability  density  function  (66)  yields  [5,  6.611  9] 


Ph(0)  -i  arc  cos(j^)  ,  1-Ph(0)  -  i  arc  cos(^) 


(69) 


23 


TR  7045 


DISTRIBUTIONS  FOR  N=3 

Use  of  [6,  10.2.17]  on  (60)  immediately  yields  probability  density 
function 


Ph(u)  - 


zum exp(u) 


'l+R)  exp^l+2R/ 


for  u  <  0 


for  u  >  0 


^  2(1+R) 

The  cumulative  and  exceedance  distribution  functions  easily  follow  as 


Ph(u)  -  2(WeXp(u) 

1  ‘  ph(u)  "  2(M)  exp(t+a!') 

This  latter  result  corroborates  (61)  and  (62). 


for  u  <  0  , 
for  u  >  0  . 


DISTRIBUTIONS  FOR  N=4 


(70) 


(71) 


The  only  closed  form  result  possible  is  obtained  by  direct  integration  of 
probability  density  function  (60)  to  get  origin  value 


n  if  /  R  \  R(1+2R)1/2 

ph(°)  =  7  [arc  cos(p£)  -  ~^2~ _ 


(72) 


This  follows  by  use  of  the  integral 
dx  e-ax  X  K.(bx)  . 

(B2-a2) 


- x-a  -x  for  B  >  -a  , 

B  (B-a) 


(73) 


which  follows  from  [5,  6.611  9]  by  applying  a/ae  to  both  sides. 
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DISTRIBUTIONS  FOR  N=5 

Use  of  [6,  10.2.17]  on  (60)  immediately  yields  probability  density 
function 


1+2R  (1+R)u  eXp(u)  for  u  <  o 


Ph(u)  -< 


4(1+R)' 


1+2R+(1+R)u 


4(1+R) 


exp 


G+2r) 


for  u  >  0 


The  cumulative  and  exceedance  distribution  functions  follow  as 


(74) 


P  (u)  =  2+3R  (1*R)u  exp(u)  for  u  <  0, 

"  4(1+R)3 

1  ~  Ph(U)  =  4~('l +R ) 3  [(1+2R)(2+R)+(1+R)u]  exp(w)  for  u  >  0  .  (75) 

This  latter  result  corroborates  (61)  and  (62).  Also,  this  example  was  used  as 
a  check  on  the  numerical  procedure  [3]  applied  directly  to  the  characteristic 
function,  which  is  used  in  the  following  section;  the  agreement  was  ten 
decimals  for  numerous  values  of  R  and  u. 
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GRAPHICAL  RESULTS  FOR  y=l,  SAMPLE  MEAN  REMOVAL 
SUMMARY  OF  PARTICULAR  CASE  CONSIDERED 


We  first  summarize  here  the  particular  case  that  will  be  considered 
quantitatively  in  this  section.  The  input  samples  are 


un  =  uu  +  us<")  +  ud(") 

vn  -  +  vs(n)  +  vd(n) 


for  1  <  n  <  N, 


where  these  Gaussian  random  variables  have  statistics 


(76) 


us(n)  -  vs(n)  =  ud(n)  =  vd(n)'  =  0  , 


u2(n)  =  S  ,  v2(n)  =  S  ,  u  (n)v  (n)  =  (S  S  )1/2  , 
sx  US  V  S  S  v  U  V 


ud(n)  =  Du,  v2(n)  =  Dy,  ud(n)vd(n)  =  0  .  (77) 

We  presume  that  the  simultaneous  signal  components  u$(n),  v  (n)  in  the  two 
channels  are  fully  correlated,  that  all  other  random  variables  are 
independent,  and  that  the  two  channel  input  signal-to-noise  ratios 


TT  =  JT 
u  v 


=  R 


(78) 


have  a  common  value  R.  More  general  situations  have  been  considered  in 
earlier  sections;  however,  only  this  special  case  will  be  numerically 
evaluated  here. 
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The  normalized  crosscorrelator  output,  with  sample  mean  removal  (Y=l),  is 


h  = 


(D  D  )112  ^ 
'  u  v7  n=l 


55  J 

n  n 


where  the  sample  ac  components 

N  N 

“n  -  un  -  I  5  um  '  *n  =  vn  ‘  ¥  ]E  vn,  • 

m=l  m=l 

The  characteristic  function  of  h  is  given  by  (57)  as 


fh(J)  =  [(1+15) (l-i|(l+2R) )] 


(79) 


(80) 


(81) 


and  depends  only  on  signal-to-noise  ratio  R  and  number  of  terms  N.  We  must 
have  N>2. 


If  R=0  and  we  evaluate  the  exceedance  distribution  function  corresponding 
to  (31),  we  then  have  the  false  alarm  probability.  But  when  R>0,  the 
exceedance  distribution  function  corresponding  to  (81)  is  the  detection 
probability.  In  the  following,  we  plot  the  detection  probability  vs.  the 
false  alarm  probability,  with  signal-to-noise  ratio  R  as  a  parameter; 
different  values  of  N  are  handled  in  separate  plots. 

OPERATING  CHARACTERISTICS  FOR  Y-1 

A  sample  program  for  evaluation  of  the  cumulative  and  exceedance 
distribution  functions  corresponding  to  characteristic  function  (81),  and 
thereby  the  detection  probability  vs.  false  alarm  probability  operating 
characteristics  of  the  crosscorrelator  with  sample  mean  removal,  is  given  in 
appendix  D.  It  is  heavily  based  on  the  technique  developed  and  explained  in 
[3]. 


28 


TR  7045 


In  figures  1-14  are  presented  the  operating  characteristics  for  the 
crosscorrelator  with  sample  mean  removal,  for  values  of 

N  =  2,  3,  4,  6,  8,  12,  16,  24,  32,  48,  64,  96,  128,  256,  (82) 

respectively.  The  case  of  N=2  was  accomplished  by  use  of  ( 67 ) — ( 69 ) ;  results 
for  N=3  relied  on  (71);  and  the  remainder  for  N>4  employed  a  numerical 
procedure  [3]  proceeding  directly  from  characteristic  function  (81)  to  the 
exceedance  distribution  function.  False  alarm  probabilities  Pp  in  the  range 
1 E— 10  to  .5  and  detection  probabilities  Pp  covering  1 E— 10  to  .999  are 
presented.  The  abscissa  and  ordinate  on  these  plots  are  according  to  a  normal 
probability  transformation,  as  explained  below.  Values  of  signal-to-noise 
ratio  R  are  taken  as  R=2n,  where  n  assumes  values  appropriate  for  each  plot 
in  order  to  cover  the  full  range  of  probabilities  of  interest. 

GAUSSIAN  APPROXIMATION 

Suppose  the  decision  variable  of  a  processor  is  Gaussian  with  mean  and 
standard  deviation  mo>  0Q  respectively  when  the  input  signal  is  absent, 
and  m^,  when  signal  is  present.  Then  for  threshold  .A,  the  false  alarm 
probability  and  detection  probability  are 


respectively,  where  <t>  and  $  are  the  normalized  Gaussian  probability  density 
function  and  cumulative  distribution  function: 

u 

tf(u)  =  (2ir )~1/2  exp(-u2/2),  $(u)  =  J dt  <4(t)  .  (84) 

If  we  let^be  the  inverse  function  to$>,  and  define 

x  =  $j(Pp) >  y  =  $r(PD)  > 


*v5ee  paje  3" I  ei  se^ 


(85) 
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then  threshold  A  can  be  eliminated  from  (83)  to  yield 


y  = 


ml  -  mo  +  V 

°1 


(86) 


Equation  (85)  corresponds  to  the  transformation  to  normal  probability 
coordinates;  thus  a  plot  of  Pp  vs  Pp  on  normal  probability  paper  is  the 
straight  line  (86)  when  the  decision  variable  is  Gaussian  under  both 
hypotheses  of  signal  absent  as  well  as  present. 


Reference  to  (59)  reveals  that,  for  our  application. 


mQ  =  0,  m1  =  (N-l)R,  a2  =  N-l,  a2  =  (N-l) (1+2R+2R2)  ,  (87) 

since  setting  signal-to-noise  ratio  R=0  corresponds  to  hypothesis  0,  signal 
absent.  Substitution  in  (86)  yields 


y  = 


(N-l)1/2  R  +  x 


(1+2R+2R  ) 


TJT 


(88) 


that  is,  if  normalized  crosscorrelator  output  h  were  Gaussian,  the  operating 
characteristics  would  be  straight  lines  dictated  by  (88).  These  straight 
lines  are  superposed  as  dashed  lines  in  figures  12-14  for  N=96,  128,  256 
respectively.  Despite  the  large  value  of  N=96  in  figure  12,  the  Gaussian 
approximation  is  not  that  accurate,  especially  for  small  false  alarm 
probabilities  and  large  detection  probabilities.  The  exact  curve  (solid)  and 
the  Gaussian  approximation  (dashed)  cross  each  other,  and  are  labelled  at  the 
crossing  with  the  corresponding  value  of  n  in  signal-to-noise  ratio  R=2n. 

For  N=256  in  figure  14,  agreement  is  better  and  the  Gaussian  approximation  is 
probably  adequate  for  larger  N.  If  not,  an  additional  term  or  two  in  an 
Edgeworth  expansion  could  be  investigated  with  the  aid  of  the  cumulants  given 
in  (58). 
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An  obvious  shortcoming  of  the  Gaussian  approximation  (88)  may  be  seen 
immediately: 


(89) 


<  1  for  any  x  . 


Reference  to  (85)  then  yields  the  interpretation 


(90) 


for  any  Pp  . 


That  is,  as  input  signal-to-noise  ratio  R  tends  to  infinity,  the  approximate 
detection  probability  saturates  at  a  value  less  than  1,  regardless  of  the 
false  alarm  probability.  Thus  the  Gaussian  approximation  must  certainly 
deteriorate  for  large  R;  the  exact  discrepancy  for  probabilities  of  practical 
interest  is  displayed  in  figures  12-14. 
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ANALYTIC  RESULTS  FOR  Y=0,  SAMPLE  MEAN  NOT  REMOVED 


In  this  section  and  the  next,  attention  will  be  confined  solely  to  the 
case  of  scale  factor  y=0.  The  characteristic  function  of  the  crosscorrelator 
output  q  is  then  given  in  (31)  and  (32),  and  the  mean  and  variance  of  q  are 
listed  in  (33).  Due  to  the  complexity  of  the  exponential  term  in 
characteristic  function  (31),  there  are  no  general  probability  density 
function  or  cumulative  distribution  function  results  for  arbitrary  N,  like 
those  given  earlier  in  (40),  (41),  and  (44)  for  y=l*  Here,  we  can  have  N>1. 

SPECIALIZATION  TO  THE  SIGNAL  AND  NOISE  MODEL 

For  the  model  presented  earlier  in  (3)-(8),  the  original  parameters  in 
(24)  take  the  form  already  given  in  (37),  but  now  with  y=0.  Specifically, 
characteristic  function  (31)  is 


(91) 


where 


The  mean  and  variance  of  crosscorrelator  output  q  follow  from  (38)  according  to 


(93) 


33 


TR  7045 


NORMALIZED  CROSSCORRELATOR  OUTPUT 

As  in  (49),  and  for  the  same  reasons,  we  define  a  normalized 
crosscorrelator  output,  relative  to  the  additive  random  noise  disturbances 
ucj(n)  and  vcj(n)  in  (3)  and  (5),  according  to 


h  = 


<W 


172  * 


(94) 


The  characteristic  function  of  h  is  available  from  (91)  and  (92): 


=  [l-i?2a+£2(B2-a2)]  exp 


-  fq(?'<w1/2)  - 

-N/2 


fi? - 211 

L  l-iJ2a+f 


a+i  5  b 


, ^2 / . 2  2  \ 
+r  (B  -a  )J 


(95) 


where 


«  =  »S(W1/2  .  S  =  [(1*RU)(1*RV)]1/2 


a  .  NruV  b  -i«[(lH-Ru)r2  *  (1+Rv)r2  -  2»s(RuRv)1,2Vv]  • 


(96) 


and  where  we  have  defined 


n  Up  V  ^U  -V 

u  =  D u  *  Rv  -  D v  ’  ru  -  1/2  ’  rv  -  1/2  • 

U  V 


(97) 


The  characteristic  function  in  (95)  depends  on  six  fundamental  parameters. 

namely  N,  p$>  ru>  rv,  ru,  rv.  The  mean  and  variance  of  h  follow  from 
(94),  (93),  and  (97): 
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“h  ■  Nf»s(RuRv)1/2  +  rurJ  ■ 

°h  •  NC(1+Ru>(1+Rv)+‘,sRuRv*(1*Ru)rv*(1*Rv)ru*2os(RuRv)1/2  Vv3  •  (981 

and  rv  are  referred  to  as  normalized  means. 

REDUCTION  TO  IDENTICAL  SIGNAL  COMPONENTS 


In  order  to  prepare  for  numerical  evaluation  of  the  operating 
Characteristics  of  the  crosscorrelator  with  y=0,  we  further  specialize  the 
signal  and  noise  model  to  the  case  where 


R  =  R  =  R,  0  =  1,  r  =  r  =r; 

u  v  ’  s  ’  u  v 


see  (56)  et  seq.  This  leads  to 


a  =  R,  8  =  1+R,  a  =  b  =  Nr  , 


via  (96).  The  characteristic  function  in  (95)  then  reduces  to 


fh(F)  -  C(l+i?)(l-ir(l+2R))]-N/2  exp 


2  T 


iJNr 

1— i V( 1+2R) J  ’ 


and  the  mean  and  variance  in  (98)  become 


uh  =  N(R+r2),  o2  =  N[l+2R+2R2+2(l+2R)r2]  . 


(99) 


(100) 


(101) 


The  characteristic  function  in  (100)  has  a  branch  point  at£=  i,  and  another 
branch  point  at  f  =  -i/(l+2R)  which  overlaps  an  essential  singularity;  this 
complicates  some  of  the  analytical  development  to  follow. 


There  are  three  fundamental  parameters  in  (100),  namely  N,  R,  r.  Since 
normalized  mean  r  appears  only  through  its  square,  we  can  presume  r>0  without 
loss  of  generality.  Furthermore,  if  r=0,  characteristic  function  (100) 
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reduces  to  (57)  if  N-l  there  is  replaced  by  N.  Thus  the  curves  for  r=0  here 
can  be  obtained  from  the  earlier  curves  for  y=l  in  figures  1-14  by  looking  at 
a  value  for  N  there  which  is  one  greater;  accordingly  we  can  confine  attention 
to  r>0  in  this  and  the  next  section. 

When  the  random  signal  components  u$(n),  v$(n)  in  model  (3)  are 
absent,  then  R=0,  and  the  exceedance  distribution  function  corresponding  to 
characteristic  function  (100)  becomes  the  false  alarm  probability.  However, 
(100)  still  depends  on  r,  meaning  that  the  false  alarm  probability  must, 
likewise.  Thus,  a  non-zero  mean  in  (3),  i.e.  r>0,  is  not  considered  a  signal 
attribute  here,  but  rather  is  a  nuisance  quantity;  it  may,  in  fact,  degrade 
the  operating  characteristics  if  not  removed. 

For  notational  convenience  in  the  following,  we  define 


(102) 


a)  =  1+2R  . 


The  magnitude  of  the  exponential  term  in  characteristic  function  (100)  then 
can  be  expressed  as 


(103) 


l+JW 


which  is  monotonical ly  decreasing  for  £ >0.  Coupled  with  the  observation 
immediately  under  (39),  it  is  seen  that  |fh(f)|  in  (100)  is  monotonical ly 
decreasing  for  all  y>0  and  any  N,R,r.  This  property  allows  for  a  convenient 
termination  procedure  in  the  numerical  transformation  [3]  of  characteristic 
function  (100).  It  should  however  be  observed  that  (103)  does  not  decrease  to 
zero,  but  saturates  at  value  exp(-Nr$),  regardless  of  how  f  increases  to 
infinity;  thus  the  eventual  decay  of  the  characteristic  function  (100)  is 
furnished  only  by  the  leading  factor. 
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ASYMPTOTIC  BEHAVIOR  OF  CUMULATIVE  AND  EXCEEDANCE  DISTRIBUTION  FUNCTIONS 

In  appendix  E,  it  is  shown  that  if  characteristic  function  (100)  is 
substituted  in  (63)  and  (64),  and  the  contours  moved  appropriately  in  the 
complex  "£-pl ane,  then  the  following  asymptotic  behaviors  obtain.  The 
cumulative  distribution  function 


2N/2<i,r)N/2]-1 

_  N(1+2R)  ,  Nr2 

wj. 


exp 


[■ 


Nr 


,21 

wj 


as  u-»  - °o. 


(104) 


We  see  again,  in  similar  fashion  to  (61)  and  (62),  that  in  order  for  the 
correction  term  in  the  second  line  of  (104)  not  to  be  too  significant,  we  must 
have  u  <  -N^.  For  reasons  elucidated  in  (61)  et  seq.,  (104)  is  not  useful 
for  large  N. 


As  checks  on  (104),  we  note  that  for  r=0  and  N=2,  (104)  reduces  precisely 
to  the  upper  line  of  (71);  this  latter  result  pertains  to  y=l>  N=3  and  is 
consistent  with  the  observation  already  made  in  the  paragraph  below  (101).  In 
addition,  if  we  let  r=0  and  N=4  in  (104),  it  reduces  to  the  upper  line  of 
(75);  this  latter  result  holds  for  y=l»  N=5  and  is  likewise  consistent. 


Also  given  in  appendix  E  are  a  variety  of  asymptotic  expansions  for  the 
exceedance  distribution  function;  the  simplest  one  is 


1  -  Ph(u) 


N-l 


2l,2  2N/2  (1+r)N/2  (Nr2) 


-1 


(1+2R) 


N+l 

~T 


N-3 

4 

u  exp 


'  fu1'2  -  N1/2rV 

■  — Fzr — 


as  u  -*  +“> ;  r>0 


(105) 
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However,  the  same  reservations  as  above,  regarding  u  large  relative  to  N^, 


are  again  in  order. 

DISTRIBUTIONS  FOR  N=1 

If  characteristic  function  (100)  with  N=1  is  substituted  in  (63),  and  if 
the  contour  is  moved  as  indicated  under  (66),  there  follows  the  exact  result 
for  the  cumulative  distribution  function 


ph<u>-f 


2  2  1/2 
0J  ( l+t2 )  ( l+w(  l+t2 ) ) 


(This  reduces  to  (67)  for  r=0,  as  it  must.)  This  integral  form  possesses  all 
the  desirable  attributes  listed  under  (67). 

If  characteristic  function  (100)  with  N=1  is  substituted  in  (64)  in  an 
attempt  to  get  the  exceedance  distribution  function,  the  analysis  becomes 
rather  difficult,  due  to  the  overlapping  essential  singularity  and  branch 
point  of  the  integrand  at  ^  =  -i/io.  This  problem  is  treated  in  detail  in 
appendix  F,  with  the  result  that  the  exceedance  distribution  function  can  be 
found  via  the  characteristic  function  approach  in  terms  of  two  integrals;  see 
( F— 21 )— (F— 23) .  However,  a  better  numerical  procedure  for  the  exceedance 
distribution  function  is  the  direct  result  derived  in  (F-33);  this  latter 
integral  is  the  one  actually  used  here  to  generate  the  operating 
characteristics  for  y=0,  N=l. 

DISTRIBUTIONS  FOR  N=2 

The  characteristic  function  is  available  from  (100): 
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where  <u  =  1+2R.  The  probability  density  functions  corresponding  to  these  two 
characteristic  functions  are  [5,  6.631  4] 

p^u)  =  exp(u)  U ( — u )  , 

P2(u)  .  i  exp(  Ip^Zu)1'2)  U(u)  ,  (108) 

where  U  is  the  unit  step  function 

U(u)  - 

The  probability  density  function  of 

+o« 

ph(u)  =  J dt  p^(t)  p^u-t)  for  all  u  .  (HO) 


for  u<0 
for  u>0 


(109) 


h  is  given  by  convolution 


Substitution  of  (108)  in  (110)  yields 


where  the  Q-function  is  defined  in  [7]  and  the  two  integrals  encountered  have 
been  evaluated  by  use  of  [5,  6.631  4]  and  [7,  (9)],  respectively.  The 
cumulative  and  exceedance  distribution  functions  of  h  readily  follow  from 
(111),  the  latter  by  means  of  [7,  (42)]: 
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ph(u)  =  i^exp(u  -t£)  for  u  1  0  » 


1  -  Ph(u)  =  Q 


k /21A 

l0,1/2  *  1“  / 


1/21 


2  r  /2(l+a>)u\ 

TC  CA*T  ■  if;  4(  T777FTT77  ’  (  S  V 


1  exp(u  -  S-)  Q/ 


1/21 


for  u  >  0.  (112) 


(As  a  check,  for  r=0,  then 
respectively,  as  they  must.) 

DISTRIBUTIONS  FOR  N=4 


0)  '  (1+0))' 

(111)  and  (112)  reduce  to  (70)  and  (71) 


The  characteristic  function  is  available  from  (100): 


fh(F)  -  (l+if )“2*(l-iJo>)-2  exp(|^j  =  ^C?)*^) 


(113) 


where  «>  =  1+2R.  The  probability  density  functions  corresponding  to  these  two 
characteristic  functions  are  [5,  6.631  4] 

P^u)  =  -u  exp(u)  U ( — u )  , 

p2(u)  =  exp(~  U+^~/)  ul/2  Il(4rui/2/a))  u(u)  •  (H4) 

The  probability  density  function  of  h  is  given  by  convolution  (110).  In 
preparation  for  that  result,  we  use  the  shorthand  notation 


5„<“>  -  &n  •  &D- 


(115) 
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where  the  Q^-f unction  is  defined  in  [8].  We  also  present  a  new  integral 
result  that  will  be  needed  in  the  sequel, 

J  dx  x  QM(b,  ax)  =  -^b2QM+i(b,ac)+2M  Qr/,(b,ac)-a2c^QM_1(b,ac)],  (116) 

which  can  be  interpreted  as  the  limit  of  [8,  (31)]  as  p-»0+.  Substitution  of 
(114)  in  (110)  yields  probability  density  function 


r 


ph(u>  - 


(l+u) 


4  exp  lu 


(1+u) 


4  expfu 


-  Y+~j  j4r2+2d)(l+a))-(l+o))2uJ  for  u  <  0 

-  •p-jpr2q3(u)+2a)(l+a))q2(u)-(l+a))2uq1(u)J  for  u  >  0 


(117) 


where  (115)  has  been  used;  the  upper  line  employed  [5,  6.631  4],  while  the 
lower  line  used  (116)  and  an  integration  by  parts  procedure  to  be  elaborated 
below  in  the  exceedance  distribution  function  evaluation. 


The  cumulative  distribution  function  for  u  <_  0  follows  readily  by 
integration  of  (117): 

P^u)  = — " — 4  exp/u  -  j^J[4r2+(l+u)(l+3a))-(l+(jl))2uJ  for  u  £  0.  (118) 

For  u  >_  0,  a  modified  approach  is  required.  Integration  of  (110)  yields 
cumulative  distribution  function 

+  «4 

Ph(u)  =  f dt  P2(t)  Pi (u-t)  .  (119) 

Probability  density  functions  and  p2  are  available  in  (114),  and 
cumulative  distribution  function  follows  as 
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Px(u)  = 


(l-u)exp(u)  for  u  £  O^j 


1  for  u  >  0 


(120) 


Substitution  of  (114)  and  (120)  in  (119)  yields,  for  u  >  0, 


ph(u)  -  £  dt  2* 


K'2' 


l(4rt1/2/„) 


dt  exp^-  — -jr-  ^  t^2  I^rt1'2/.)  (1-u+t)  exp(u-t)  = 
u 


+O0 

=  l-q2(u)  +2^  exp  |u-  — j"dt  t1/2  exp^-  t  I^rt1^/®)  (1-u+t),  (121) 


via  [8,  (22)].  We  now  integrate  by  parts,  letting  U(t)=l-u+t,  and  the 
remainder  dV(t).  Then  using  [8,  (22)]  again,  we  find 


V(t)='^eXP(™q2(t) 

Combining  these  results,  (121)  and  (116)  yield 


(122) 


=  ^ 


1  -  Ph(u)  »  g2(u)  -  exp^u  -  q2(u)  -  J"dt  q2(t)  . 

(u) - ^  exp ^  -  yy[4i-2qj(u)'-(l'-»)(l+3»)q2(u)-(l+u)2uq1(u] 


for  u  >  0.  (123) 
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The  final  results  for  N=4  are  given  by  (115),  (118),  and  (123).  This  case  was 
used  as  a  numerical  check  on  the  computational  approach  [3]  proceeding 
directly  from  characteristic  function  (100)  to  the  exceedance  distribution 
function,  with  excellent  agreement  for  numerous  values  of  R,  r,  and  u. 
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GRAPHICAL  RESULTS  FOR  Y=0,  SAMPLE  MEAN  NOT  REMOVED 

SUMMARY  OF  PARTICULAR  CASE  CONSIDERED 

The  situation  of  interest  has  already  been  summarized  in  (76)— (78) ;  in 
addition,  we  have  a  common  value  for  the  normalized  means. 


0l/2  *  „l/2  *  r  ' 


(124) 


and  the  normalized  crosscorrelator  output  is  not  (79)— (80) ,  but  rather  is,  for 

Y-0, 

N 

(125) 


h  = - - — i-pr  >  u  v  . 

(r\  n  \l n  n 


<W  n=l 


The  characteristic  function  is  given  by  (100): 


v?> 


[(l+if)(l-ifa>)rN/2 


(126) 


where  u>  =  1+2R.  When  the  signal  is  absent,  then  R=0;  however  the  false  alarm 
probability  corresponding  to  this  characteristic  function  still  depends  on  r. 
Thus  since  the  sample  mean  has  not  been  removed,  the  operating  characteristics 
will  also  depend  on  r.  Since  results  for  r=0  can  be  found  from  an  earlier 
section,  we  only  consider  r>0  here. 


OPERATING  CHARACTERISTICS  FOR  Y=0 

A  sample  program  for  evaluation  of  the  cumulative  and  exceedance 
distribution  functions  corresponding  to  characteristic  function  (126),  and 
thereby  the  detection  probability  vs.  false  alarm  probability  operating 
characteristics  of  the  crosscorrelator  without  sample  mean  removal,  is  given 
in  appendix  G.  In  figures  15-35  are  presented  the  operating  characteristics 
for  values  of 

.Te*  se< 
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N  =  1,  2,  3,  4,  8,  16,  32,  64,  128,  256,  (127) 

and  for  various  values  of  r.  The  case  of  N=1  was  accomplished  by  use  of  (106) 
and  (F-33);  results  for  N=2  employed  (112);  and  the  remainder  for  N>3  employed 

a  numerical  procedure  [3]  proceeding  directly  from  characteristic  function 

(126)  to  the  exceedance  distribution  function.  False  alarm  probabilities  Pp 
in  the  range  IE— 10  to  .5  and  detection  probabilities  Pp  covering  IE— 10  to 
.999  are  presented.  The  abscissa  and  ordinate  on  these  plots  employ  a  normal 
probability  transformation,  as  explained  earlier  in  (84)-(86). 

Values  of  signal-to-noise  ratio  R  are  taken  as  R=2n,  where  n  assumes 

values  appropriate  for  each  plot  in  order  to  cover  the  full  range  of 
probabilities  of  interest.  VaTues  of  normalized  mean  r  in  (124)  have  been 
taken  as  r=l  and  2,  with  the  exception  of  figure  19  where  one  example  for  r=4 
was  added. 

Without  exception,  increasing  r  from  zero  degrades  the  operating 
characteristics  of  the  crosscorrelator.  For  example,  figures  17,  18,  19  give 
a  succession  of  operating  characteristics  for  r  =  1,  2,  4  respectively,  and 
for  common  values  of  signal-to-noise  ratio  R.  (in  order  to  determine  the 
operating  characteristics  for  r=0  here,  we  can  look  at  the  earlier  results  in 
figures  1-14,  but  for  a  value  of  N  which  is  one  greater  there.)  Thus,  not 
removing  the  sample  mean  from  the  crosscorrelator  output  requires  a  larger 
threshold  setting  for  a  specified  false  alarm  probability  and  thereby  lowers 
the  detection  probability  and  degrades  performance. 

GAUSSIAN  APPROXIMATION 

If  the  crosscorrelator  output  is  Gaussian,  for  both  signal  absent  as  well 
as  present,  the  earlier  derivations  in  (83)-(86)  pertain.  Now  reference  to 
(101)  yields  statistics 


(Hq  =  Nr 


l,  =  N(R+r  ) 
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a2  =  N[l+2R+2R2+2(l+2R)r2], 
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since  setting  signal-to-noise  ratio  R=0  corresponds  to  hypothesis  0,  signal 
absent.  Substitution  in  (86)  yields  the  normal  probability  approximation 


y  = 


Nl/2  R+(l+2r2) 

[i+2R+2R2+2(l+2R) 


1/2 


x 

172  ’ 


(129) 


These  straight  lines  are  superposed  as  dashed  lines  in  figures  32-35  for  N=128 
and  256.  The  Gaussian  approximation  is  moderately  good  for  large  N  such  as 
256,  and  in  fact  crosses  the  exact  curves  (solid)  at  a  point  which  is  labeled 
with  the  corresponding  value  of  n  in  signal-to-noise  ratio  R=2n. 


An  obvious  shortcoming  of  the  Gaussian  approximation  (129)  is  apparent: 


11m  y  = 
R-*°6 


Reference  to  (85)  then  yields  the  interpretation 


r 


for  any  x,  r 


(130) 


(131) 


That  is,  as  input  signal-to-noise  ratio  R  tends  to  infinity,  the  approximate 
detection  probability  saturates  at  a  value  less  than  1,  regardless  of  the 
false  alarm  probability  and  normalized  mean  r.  Thus  the  Gaussian 
approximation  must  certainly  be  inaccurate  for  large  R;  the  exact  discrepancy 
for  probabilities  of  practical  interest  is  displayed  in  figures  32-35. 
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SUMMARY 

A  closed  form  expression  for  the  characteristic  function  of  the  output  of 
a  crosscorrelator,  with  or  without  sample  mean  removal,  has  been  derived  in 
(23 ) — ( 24 )  for  general  values  of:  the  number  of  terms  added  to  yield  the 
correlator  output,  the  means  and  variances  in  each  of  the  two  input  channels, 
the  degree  of  correlation  between  the  two  channels,  and  the  scale  factor 
employed  in  the  sample  mean  removal.  A  program  for  the  evaluation  of  the 
cumulative  and  exceedance  distribution  functions  of  this  general  case  has  also 
been  presented.  These  results  can  furnish  the  basis  of  a  study  of  the  error 
probabilities  of  a  correlator  required  to  decide  between  alternative 
hypotheses  on  the  input  statistics  [1];  this  problem  will  in  fact  be  the 
subject  of  a  future  technical  report  by  this  author. 

The  general  results  were  first  specialized  to  a  signal  and  noise  model, 
and  then  to  the  two  distinct  cases  of  sample  mean  removal  (y=l)  and  no  sample 
mean  removal  (y=0).  Plots  of  the  operating  characteristics  for  numerous 
values  of  N  and  signal-to-noise  ratio  R  were  then  displayed  for  a  wide  range 
of  detection  probability  vs  false  alarm  probability.  Some  new  analytic 
results  for  cumulative  and  exceedance  distribution  functions,  especially  for 
small  N,  were  derived  and  used  as  checks  on  the  general  numerical  procedure. 
Comparisons  with  a  Gaussian  approximation  indicated  quantitatively  when  that 
simplification  is  valid.  Asymptotic  results  derived  were  useful  for  small  N, 
but  not  for  large  N  except  in  the  region  of  probabilities  too  small  to  be  of 
practical  importance. 
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IE- 10  1 E-8  1 E-B  1 E-4  .001  .01  .1  .5 


Figure  1.  Operating  Characteristics  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  2 
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IE- 10  IE-8  IE-6  IE-4  .001  .01  .1  .5 


Figure  2.  Operating  Characteristics  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  3 
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IE- 10  IE-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  3.  Operating  Characteristics  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  4 
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Figure  4.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  6 
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IE- 10  IE-8  1 E-B  IE-4  .001  .01  .1  .5 


P 

F 

Figure  5.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  8 
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IE-10  1 E-8  IE-6  1 E-4  .001  .01  .1  .5 

P  * 

F 

Figure  6.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  12 
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IE-10  IE-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  7.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  1 6 
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IE- 10  IE-8  1 E-6  1 E-4  .001  .01  .1  .5 

P 

F 

Figure  8.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  24 
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IE- 10  IE-8  IE-6  1 E-4  .001  .01  .1  .5 


P 

F 

Figure  9.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  32 
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IE- 10  1 E-8  1 E-G  IE-4  .001  .01  .1  .5 


Figure  10.  Operating  Ch ar act e r i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  48 
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IE- 10  1 E-8  IE-6  IE-4  .001  .01  .1  .5 


Figure  11.  Operating  Characteristics  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  64 
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IE-10  1 E— 8  IE-6  1 E-4  .001  .01  .1  .5 


Figure  12.  Operating  Character i st i cs  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  -  96 
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IE- 10  IE-8  IE-6  IE-4  .001  .01  .1  .5 


? 


Figure  13.  Operating  Char acter i st i cs  for 
Correlator  with  Sample  Mean  Removal,  N 


Cross- 
=  128 
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IE- 10  1 E-8  1 E-B  1 E-4  .001  .01  .1  .5 


Figure  14.  Operating  Characteristics  for  Cross- 
Correlator  with  Sample  Mean  Removal,  N  =  25G 
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IE- 10  IE-8  1E-G  1 E-4  .001  .01  .1  .5 


Figure  15.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N= 1 ,  r=l 
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IE- 10  1 E-8  1 E-B  IE-4  .001  .01  .1  .5 

D 

V 

Figure  16.  Operating  Char acter i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=l,  r=2 
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IE- 10  IE-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  17.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=2 ,  r=l 
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IE- 10  IE-8  1E-B  IE-4  .001  .01  .1  .5 


Figure  18.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=2  ,  r=2 
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IE- 10  IE-8  1 E-G  1 E-4  .001  .01  .1  .5 


Figure  19.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=2 ,  r=4 
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IE-10  IE-8  1 E-G  IE-4  .001  .01  .1  .5 


Figure  20.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=3 ,  r=l 
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IE- 10  IE-8  1 E-B  IE-4  .001  .01  .1  .5 


Figure  21.  Operating  Ch ar acte r i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=3 ,  r=2 
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Figure  22.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=4 ,  r=l 
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IE-10  1 E-8  1 E-B  IE-4  .001  .01  .1  .5 


Figure  23.  Operating  Ch ar act e r i st i cs  -for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=4 ,  r=2 
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IE-10  IE-8  1E-B  1 E— 4  .001  .01  .1  .5 


Figure  24.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=8 ,  r=l 
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IE- 10  IE-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  25.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=8 ,  r=2 
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IE- 10  IE-8  IE-6  1 E-4  .001  .01  .1  .5 


Figure  26.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=  1  6  ,  r=l 
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IE- 10  IE-8  1 E-B  IE-4  .001  .01  .1  .5 


Figure  27.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=  1  6  ,  r=2 
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IE- 10  1 E-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  28.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=32 ,  r=l 
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IE-10  IE-8  1 E-G  IE-4  .001  .01  .1  .5 


Figure  29.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=32 ,  r=2 
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Figure  30.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal ,  N=64 ,  r=l 
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IE- 10  IE-8  1E-B  IE-4  .001  .01  .1  .5 


Figure  31.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=64 ,  r=2 
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IE- 10  IE-8  IE-6  IE-4  .001  .01  .1  .5 


Figure  32.  Operating  Character i st i cs  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=128,  r=l 
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IE-10  IE-8  IE-6  IE-4  .001  .01  .1  .5 


Figure  33.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=128,  r=2 
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IE- 10  IE-8  1 E-6  IE-4  .001  .01  .1  .5 


Figure  34.  Operating  Characteristics  for  Cross- 
Correlator  without  Sample  Mean  Removal,  N=256 ,  r  «=  1 
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Figure  35.  Operating  Character i st i cs 
Correlator  without  Sample  Mean  Removal, 


for  Cross— 
N=256 ,  r =2 

■J^sve-vse/  bfflnk 
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APPENDIX  A.  CORRELATOR  OUTPUT  INDEPENDENCE  OF  MEANS 

If  we  let  scale  factors  a=l  but  s^l  in  (11),  we  still  get  y=1  from  (13). 
This  means  that  correlator  output  q  in  (12)  is  independent  of  uu,  „  g. 

To  see  this  directly,  let 

un  =  uu  +  ^n’  vn  =  uv  +  zn  ’  (A-l) 


where  means 


yn  -  zn  .  0  .  (A-2) 

Then  (11)  yields,  with  <*«!, 


N 


m=l 


which  is  obviously  independent  of  the  actual  value  of  mean  u  Also,  (11) 
yields 


N 

V'ril2.*  “v*1-6’  ■  <A-4) 

m=l 

which  still  depends  on  uv  and  6.  When  (A-4)  is  substituted  in  (12),  we  get 
correlator  output 


=  5 

n=l 


■n  "  I 


m 


»,(!-«) 


m=l 


(A-5) 


But  since 


N 

2E  Un  =  0  (A-6) 

n=l  n 


A-l 
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from  (A— 3 ) ,  ( A— 5 )  reduces  to 

N  N  /  N  \ 

q  =  i>u  z  =  S  (y  -  i  y  )  z  = 

n=l  n=l V  m=l  / 

N  N  N 

=  Vn  ]>2n  <A'7> 

n=l  m=l  n=l 

in  terms  of  the  ac  components  defined  in  (A-l)  and  (A-2).  Correlator  output 
( A— 7 )  is  obviously  independent  of  means  uu,  „v  and  scale  factor  b  in  (11). 


A-2. 
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APPENDIX  B.  A  USEFUL  INTEGRAL  OF  EXPONENTIALS  OF  MATRIX  FORMS 


For  symmetric  KxK  matrix  M,  with  det  M>0,  the  following  K-fold  integral 
is  well  known  (see  for  example,  [9,  section  8-3]): 


I 


dX  exp[- 


2-XTMX+NTX]  = 


( B— 1 ) 


Here  X  and  N  are  Kxl  column  matrices.  We  wish  to  extend  this  result  to  the 
case  of  double  integral 


I  =  ( 2 tt )  K  Jj  dU  dV  exp  jj-  UTAU  -  |  VTBV  + 


UTCV 


DTU 


etv 


( B— 2 ) 


where  A  and  B  are  symmetric  without  loss  of  generality,  and  the  integral 
converges;  here,  matrices  A,  B,  C  are  KxK  while  U,  V,  D,  E  are  Kxl.  Notice 
that  if  we  had  the  apparently  more-general  term 


U^V  +  VTC?U  =  UT(C1  +  cJ)V  ,  (B-3) 

we  would  simply  let  C  =  +  C^,  and  thereby  immediately  have  form  (B-2). 

To  accomplish  the  evaluation  of  I  in  (B-2),  identify  M  =  B ,  N 1  =  U^C  +  E^", 
X  =  V  in  (B-l),  and  thereby  evaluate  the  V-integral  in  (B-2),  with  result 


eWu 


E>1 


(B-4) 


Substituting  (B-4)  in  (B-2),  regrouping,  and  using  the  symmetry  of  B  (and 
therefore  B--*-),  there  follows 


I  = 


(2ir) 


— K  /  2 


(det  B) 


1/2 


dU  exp 


["  2 


ut(a-cb-1ct)u 


(DT+ETB_1CT)U 


i  etb_1e] 


(6—5) 


B-l 


TR  7045 


Now  reemploying  (B-l)  with  identifications  M  =  A  -  CB^C^,  N  =  D  +  CB  ^E, 
X  =  U,  we  get  a  closed  form  result  for  ( B— 5 ) : 


I  =  jdet(AB  -  CB-Vb)]  ^  exp|^  ETB  XE  + 

+  ^( D+CB_1E ) T ( A-CB-1CT ) _1 ( D+CB-1E )J  .  (B-6) 

This  is  the  desired  general  result  for  integral  (B-2). 


As  an  aside,  there  is  probably  a  more  symmetric  closed  form  result  than 
(B-6),  since  if  we  represent  (B-2)  by  I(A,  B,  C,  D,  E),  we  quickly  see,  by 
interchange  of  dummy  variables  U  and  V,  that 

I(A,  B,  C,  D,  E)  =  I ( B ,  A,  CT,  E,  D).  (B-7) 

However,  we  have  not  discovered  the  symmetric  form  of  (B-6).  The  present  form 
follows  as  a  result  of  the  sequential  integration  of  (B-2),  first  on  V,  then 
on  U. 


B-2 
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APPENDIX  C.  PROGRAM  FOR  CUMULATIVE  AND  EXCEEDANCE  DISTRIBUTION 
FUNCTIONS  VIA  CHARACTERISTIC  FUNCTION  (23)-(24). 

The  numerical  procedure  employed  in  appendices  C,  D,  and  G  here  is 
heavily  based  on  [3].  The  choices  of  L,  A,  b  in  lines  90  to  110  to  control 
truncation  error  and  aliasing  are  also  made  according  to  the  method  of  [3]. 
The  parameters  in  (24)  are  evaluated  once  in  lines  210-260  so  as  to  minimize 
computation  time.  The  FFT  subroutine  used  in  lines  1030  et  seq  is  listed  in 
[3,  pp.  B— 11  -  B— 12] ,  and  employs  a  zero-subscripted  array.  A  sample  plot  of 
the  cumulative  and  exceedance  distribution  functions  follows  the  program. 


10  ! 

CROSS-CORRELATOR  WITH  SAMPLE  MEAN  REMOVAL;  NUSC  TR  7 

045 

20 

N  =  5  ! 

N  u  m  b  e  r  o  f  t  er  m  s  s  u  m  m  e  d  t  o  y 

i  e  1  d  o  u  t  p  u  t 

30 

Gamma =.5  1 

Scale  factor  in  sample  mean 

removal 

40 

H  u = -  •  4  ! 

U  channel  mean 

50 

Mu*. 3  1 

V  channel  mean 

60 

S  u  = .  7  1 

U  c  hi  an  n  e  1  s  t  an  d  ar  d  d  e  v  i  at  i  o  n 

70 

3  = ,  9  i 

V  c h a n n  e 1  stand a r d  d e v i a t i o n 

8  0 

R  h  o  = .  6  ! 

C  o  r  r  e  1  at  i  o  n  c  o  e  f  f  i  c  i  e  n  t 

90 

L  =  6  0  ! 

L  i  rn  1 1  o  n  i  n  t  e g r  a  1  o  f  c  h  ar  . 

f  L4  n  c  t  i  o  n 

1  0  0 

Delta*. 12  1 

S  am  p  1  i  n  g  i  n  c  r  e  m  e  n  t  o  n  c  h  ar  . 

f  L4  n  c  t  i  o  n 

110 

Bs=. 25*<2*PI/Be 1 ta>  ! 

Shift  b,  as  fraction  of  al  i 

as  interval 

120 

Mf=2A8  ! 

S  l  ze  of  FFT 

1  3  0 

PRINTER  IS  0 

148 

PRINT  "L  = 11  ;  L  ,  "Del  ta  =  11 

; Del ta, "b  =  " ;Bs,HMf  =  "  ;  Mf 

150 

RED  I M  X  <  0  :  Mf-  1  >  ,  V  0 :  Mf  - 

1  > 

160 

D I M  X < 0 : 1 G2  3 > , Y 1 0 : 1 823 > 

170 

Su2=Su*Su 

!  C  a  1  c  u  1  a  t  e 

1 8  0 

Sv2=Sv*Sv 

1  parameters 

1 9  0 

T 1  =  1  —  G  arii  m  a 

2  0  0 

T2=T 1 *T 1 

210 

E  1  =  2  *  F:  hi  o  *  S  u  S  v 

2  2  0 

E2=S u 2  * S v 2* ( 1 - R  h o * R h o > 

230 

F 1 =E 1 *T 1 

240 

F2=E2*T2 

250 

G 1 =H*Mu*Mu*T 1 

2  60 

G  2  = .  5  *  N  *  <  S  u  2  *  M  v  *  M  u  +  S  v  2  *  M  u  *  M  u  -  E 1  *  M  u  *  M  v  >  *  T  2 

270 

N  1  = .  5  *  <  N  - 1  > 

2  S  0 

M  u  q = G 1  +  (  N  -  G  am  rn  a  )  #  R  h  o  *  S  u 

*Sv  !  Me* an  of  random 

v  ar  i  at  1  e  q 

290 

M  u  y  =  f't  u  q  +  E  s 

!  M  e  an  o  f  s  h  i  f  t  e 

d  variable  y 

300 

X<8)=0 

310 

Y < 0 )  = .  5  *  B  e  1  ta* M u y 

320 

FOR  Ms  =  1  TO  I N T < L s D e 1 t a ) 

3  3  0 

X  i  -De  1 1  a* Ns 

!  Argument  xi  of 

char,  function 

340 

X 2  =  X i *Xi 

!  C  a  1  c  u  1  at  i  o  n 

350 

T  1  =  -  X  i  *  F  1 

!  of 

360 

T2=1+X2*F2 

!  c harac t  er i st i c 

370 

CALL  In  v  <  - X 2 *  G2 ,  X  i  * G  l  ,  T 

2,T1,h,E>  1  function 

330 

C ALL  Lo g < 1 +X2 * E 2 , -X i *E 1 

,  C ,  D )  !  f’yCxi  > 

390 

CAL  L  L  tfg  <  T  2.,  T 1  ,  E ,  F  > 

400 

CALL  Exp (  R-N 1 *C- i 5*E , E- 

N1*D-. 5*F+Xi *Es, Fyr, Fy i > 

C-l 


I 


Co  1  1  dps i n g 


TR  7045 

410  Ms  =  Ns  MOD  Mf 

420  XCMs>=X<Ms)+F yr  Ns 

4  ■';!  0  Y  c.  M  s —  Y  *•.  M  s  +  F  y  i  •••  M  s 

440  NEXT  Ns 

450  CALL  Ff  MOzcMf  5  X«:>>  ,  Y  <  +  >:■ 

460  PLOTTER  IS  "GRAPHICS1' 

4  7  0  G  R  H  P  H  I  C  S 

4 8 0  SCALE  O, Mf , -14, O 

490  LINE  TYPE  3 

500  GRID  Mf/8,1 

510  PEHUP 

520  LINE  TYPE  1 

530  B=Bs*Mf*De1 1 a-' < 2*P  I  ) 

540  MOVE  B , 0 

550  DRAW  B , - 1 4 

560  PENUP 

570  FOR  K s  =  0  TO  Mf-1 

5  8  6  T  =  Y  <  K  5-  >  /  P  I  -  K  s  /  M  f 

590  X  < K s  >  = . 5 - T 

600  Y  < K s >  =  P r = . 5  +  T 

610  IF  P  r  >  = 1 E - 1 2  THEN  Y=LGT  <Pr  > 

620  IF  P  r  <  =  -  1  E  -  1 2  THEN  Y=-24-LGT < -Pr > 

630  IF  ABSCPrXlE-12  THEN  Y  =  -12 

640  PLOT  K s , Y 

650  NEXT  Ks 

660  PENUP 

670  PR  I  NT  Y  <  0 )  ;  Y  ■;  1  >  ;  Y  <  Mf -2  >  ;  Y  (.  Mf-1  > 

6  8  O  F  0  R  K  s  =  0  T  0  M  f  -  1 

690  Pr=X<Ks  > 

700  IF  P  r  >  =  1  E  - 1  2  THEN  Y  =  L  G  T  <  P r > 

710  IF  Fr<=-1E-12  THEN  Y=-24-LGT < -Pr > 

720  IF  RBSCPrX  IE- 12  THEN  Y  =  -12 

730  PLOT  Ks , Y 

740  NEXT  Ks 

750  PENUP 

760  PAUSE 

7  70  DU M P  1 3 R A P H I C S 

780  PRINT  L I N < 5 > 

790  PRINTER  IS  16 
S0O  END 

8  1 0  ! 

820  SUB  Di XI  ,  Y1  ,  X2,  Y2,  A,  B> 

830  T  =  X  2 *  X  2  +  Y  2  *  Y  2 

340  A=<X1*X2+Y1*Y2XT 

850  B  =  £  v  1  *  — X 1  &  V  2  )  *'  T 

360  SUBEND 

870  ! 

8  8  0  S  U  E  L  o  q  <  X  ,  Y ,  A ,  B  > 

890  R= . 5  *  L  0  G  C X *  X  +  Y  *■ Y  > 

900  IF  XO0  THEN  930 
910  B  = . 5 *P I *SGM < Y  > 

9 2 0  G 0 TO  95 O 

930  B  =  ATN(  YVX> 

940  IF  X<0  THEN  B  =  B  +  P I  * <1 -2* < Y< O > > 

950  SUBEND 

960  ! 


0  s  u  b  s  c  r  i  p  t  F  F  T 


Or i g i n  for  random  oar i ab  1 e*  q 


C u m u  1  at  i  o e*  p r  o  b  ab  i  1  i  t  y  i  n  X *  > 
E  x  c  e  d  a  n  c  e  p  r  o  b  a  b  i  1  i  t  y  i  n  Y  <  *  > 


PRINCIPAL  LOG  <  Z ) 


C-2 
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970  SUB  Ex p  •'  X ,  V ,  ft ,  B  )  1  E X P  (  2  > 

930  T  =  EXPOO 

990  fl=T*COS<Y) 

1000  B  =  T  *  S  I N  <  V  > 

1010  SUBENB 
1020  ! 

1030  SUB  Ff  1 1 0z  !  M  <=  2-"10  =  1024,  H=2A  INTEGER  0  SUBSCRIPT 


C-3 
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Figure  C— 1 .  Cumulative  and  Exceedance 
Distribution  Functions 


C-4 


APPENDIX  D.  PROGRAM  FOR  EVALUATION  OF  OPERATING 
CHARACTERISTICS  FOR  Y=1 
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The  comments  in  appendix  C  are  relevant  here  also.  The  characteristic 
function  used  as  the  starting  point  is  given  by  (81).  Sampling  increment  A 
employed  on  the  characteristic  function  can  be  coarse  for  small  signal-to-noise 
ratio  R,  but  must  be  finer  for  larger  R.  The  quantity  DeltaO,  aq>  in  line 
30  is  that  used  for  R=0;  all  other  A  values  are  sub-multiples,  as  indicated  by 
lines  110-130.  A  table  follows. 


N 

4 

6 

8 

12 

16 

24 

32 

48 

64 

96 

128 

256 

.10 

.09 

.08 

.06 

.05 

.05 

.05 

.05 

.04 

.03 

.03 

.02 

Table 

D-l. 

Values 

Of  A 

O 

for 

Y=1 

Let  A^  denote  the  sampling  increment  employed  for  a  particular  value  of 
signal-to-noise  ratio  R^>0.  The  cumulative  and  exceedance  distribution 
function  values  are  available  at  spacing  s  =  2tt/(M^a)  in  general,  where  Mf 
is  the  FFT  size.  If,  for  example,  Aj  =  aq/2,  then  sx  =  2sQ,  meaning 
that  probability  values  occur  twice  as  coarsely  for  R^  as  for  R=0.  Then  in 
order  to  plot  detection  probability  PD  Vs  false  alarm  probability  Pp 
without  interpolating  points,  it  is  necessary  to  skip  every  other  Pp  point 
available,  and  only  plot  those  Pp,  p^  pairs  corresponding  to  the  same 
threshold.  More  generally,  if  Aj  =  Aq/K,  where  K  is  an  integer,  then 
Sj_  =  Ksq,  and  we  plot  only  every  K-th  point  of  the  available  Pp  values. 

Here  we  have  chosen  K  to  be  a  power  of  2,  for  the  purpose  of  ease  of  plotting. 

We  choose  bias  (shift)  b  in  line  •SO  in  order  to  give  a  random  variable  y 
which  is  virtually  always  positive  for  R=0;  see  [31.  We  then  keep  b  fixed  as 
R  increases,  which  makes  the  probability  of  y>0  even  greater.  This  feature  of 
choosing  the  same  b  for  all  R^O  enables  an  easy  comparison  of  Pp  and  Pp, 
since  common  threshold  values  are  then  conveniently  realized. 
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It  was  observed  under  (39)  that  the  characteristic  functions  (39)  or  (81) 
have  monotonically  decreasing  magnitudes  for  all  £>0.  This  makes  the  choice 
of  L,  the  truncation  value  on  the  characteristic  function  integral  in  (63)  or 
(64),  rather  simple;  all  we  need  to  do  is  monitor  |f^(f)J  of  (81)  until  it 
decreases  below  a  tolerance,  here  taken  as  IE— 12 .  There  is  no  tri al-and-error 
procedure  required  as  in  [3]  to  guarantee  negligible  truncation  error. 


Subroutines  Exp  and  Log  have  already  been  listed  in  appendix  C,  and  so 
are  not  relisted  here. 

10  !  G H f'l M ft  =  1  SAMPLE  MEAN  REMOVAL 


*2  0 

Me  =8  ! 

N  ,  N  u  m  b  e  r  o  f  t,  e-  r  rn  s  s  u  rn  rn  e  d 

30 

Delta©-- 08  1 

I n i t i a 1  delta 

40 

Bs  =  P I /Be  1 t  a0  ! 

B  i  as  b 

50 

f'l  f  =  2  *‘s  1 0  i 

S i ze  of  FFT 

60 

OUTPUT  0; "GfiMMfi  =  1";"  N 

=  " ;  Nc 

70 

OUTPUT 

80 

BRTH  -2, -1,0,. 5 ,1,1. 5, 2, 2-5, 

3, 3. 5,  4, 5 

90 

RERB  Ns<*>  ! 

SNR  P  =  2 A n 

100 

OUTPUT  0;Ns<*>; 

1  1  0 

DRTR  1,2,2,2,4,4,4,8,8,16,16 

, 32, 64 

120 

R  E  R  D  I  d  e  1  t  a  (  *  ) 

130 

MRT  Belt  a® < Belt  a©  ) / 1 de 1 1 a 

1  40 

OUTPUT  O;  He  1 1  a*-  *  > ; 

150 

DRTR  IE-10,  IE-9,  IE-8,  IE-7,  IE 

-6,  1 E-5 ,  1  E-4 ,  . 001 ,  . 01 ,  .  1 ,  . 5 

160 

READ  Sc<*> 

170 

B I M  N  s  <  1 :  12),  I  d  e  1 1  a  <  © :  1 2  > ,  D  e  1 1  a  <  @ :  1 2  >  ,  s  c  <  1 :  14  > 

180 

B  I M  X  <  0:8191  > ,  V G  :  8 191) 

190 

FOR  1=1  TO  14 

200 

S'  c  <  I  >  =  F  HI  n 1  p  h  i  <  3  c  <  I  >  > 

210 

NEXT  I 

220 

S  =  8  c  (  1  > 

230 

B=Sc < 14) 

2  4  0 

Seal  e  -  <  B  -  8 )  /  03-3  ) 

250 

X  1=30 

260 

X  2  =  1  7  U 

27© 

VI  =35 

280 

Y2=Y1+<X2-X1 )*Scal e 

290 

PLOTTER  IS  "  9872R " 

300 

LIMIT  X 1 , X  2 , Y 1 , Y 2 

310 

OUTPUT  705; "VS 3" 

320 

SCRLE  S , 0 , S , B 

330 

FOR  1=1  TO  14 

340 

MOVE  S , S c  < I > 

350 

B  R  fl  N  0 ,  S  c  (.  I ) 

360 

NEXT  I 

370 

FOR  1=1  TO  11 

380 

390 

MOVE  Sc (  I)  , S 

BRflH  Sc <  I  > , B 

4  0  0 

NEXT  I 

4  1  0 

MOVE  S,S 

420 

BRAN  0,O 

4  3  0 

FEMUR 
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440  M 1 = M  f - 1 
45S  M 2  =  M f  2 

4  6  0  F  U  R  I  n =  y  TO  Ik! 

470  IF  I n > 0  T H E r \  5 0 0 
480  Rc  =0 

490  GOTO  510 

50O  R c  =  2  W s  (  I  n )  ! 

510  OUTPUT  0;"R  =  " ;Rc,"  Delta 
520  ASSIGN  #1  TO  "ABSCIS"  ! 

530  Belt a=De  1 1 a <  I n >  ! 

54  0  R  2  =  R  c *  2 

550  R2 1 =  R2  +  1 

560  N 1 2  =  <  N  c - 1 > s 2 

5  7  O  M  u  x  =  <  N  c  - 1  >  *  R  c  ! 

5  8  0  M  u  y  =  M  *.4  x  +  Bs  ! 

590  REDIM  X<0:M1>,Y<0:M1> 

600  MAT  X  =  2ER 

610  MAT  Y  =  Z E F: 

6  2  0  X  (  0  >  =  0 

6  3  0  Y  <  0  >  =  .  5  *  D  e  1  t  a*  M  u  y 

640  Ls=0 
650  L  s  -  L  s  +  1 

660  Xi=De1t.  a*Ls  1 

670  CALL  Log< 1+Xi *Xi *R21 , - X i * R 2 , 

680  CALL  Exp<-N12*Ai fXi*Bs-N12*B 

690  M s  =  L  s  MOD  Mf 

700  A  r  =  F  y  r / L  s 

710  A  i  =  F  y  i  L  s 

720  X  < M s  >  =  X  <  M s  )  +  A r 

730  V  (  M s  >  =  Y  ■'  Ms  >  +A  i 

740  M ag s  q  =  A r  *  A  r  +  A i ^A  l 

750  IF  Mags q  > 1 E-24  THEN  650 

760  OUTPUT  0 ; " X i  = " ; X i ; "  .  Mag 

770  CALL  Ff  1 13zU1f ,  X  <  *  >  ,  Y<*>  > 

780  FOR  Ms-0  TO  Ml 
790  T  =  Y ( Ms )  PI-Ms/Mf 

8  0  0  X  <  M  s  >  = .  5  -  T  ! 

810  Y  <  M  s  >  =  . 5  +  T  ' 

820  NEXT  Ms 

830  OUTPUT  0 ; Y < 0 > ; Y ( 1 > ; Y < M 1 - 1 ) ; Y 
840  PLOTTER  IS  "GRAPHICS'* 

850  GRAPHIC S 

860  SCALE  0 , M f , - 1 4 , 0 
870  LINE  TYPE  3 

880  GRID  Mf/8,1 

898  PENUP 

900  LINE  TYPE  1 

910  FOR  M s  =  0  TO  Ml 

920  P  r  =  Y  < M s > 

930  IF  Pr>=lE-12  THEN  Y  =  LGT <Pr ) 
940  IF  P r <  =  - 1 E - 1 2  THEN  Y=-24-LGT 
950  IF  ABSCPrX  IE- 12  THEN  Y  =  -12 

960  PLOT  M s  ,  Y 

970  NEXT  Ms 

980  PENUP 


S  N  R  P  =  2  ^  n 
=" ; Del ta(In) 

T  e  rn  p  o  r  a  r  y  s  t  o  r*  a  g  e  f  o  r 
f  a  1  s  e  a  1  a r  rn  fj  r  o  b  ab  i  1  i  t.  y 


M  e  a n  o  f  r  a n  d  o  rn  v  a r  i  ab  1  e  h 
t'l  e  an  o  f  s  h  i  f  t  e  d  v  ar*  i  ab  1  e  y 


A  r  g  u  m  e  n  t  x  i  o  f  c  h  a  r  .  f  r*i . 

A  i  ,  B  i  >  !  Calc  lj  1  at  i  on 

i  ,  F  y  r ,  F  y  i  *>  !  o  f 

1  characteristic 
!  f  u  n  c  t  i  o  n 
!  fy<xi> 


=  " ; SQR  < Mags q > 


C u m  '.4 1  at.  i  e  d  i  s  t  r  i  b u t  i  o n  f  u n  c  1 1  o  n 
E x c  eed  a n c  e  d  i  s t  r  i  b »j t  i  o n  f  u r*i c  t  i  o n 

<  M  1  > 


<  -Pr ) 


D-3 
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99U  FOR  M s  =  0  TO  Ml 
1 0  0  0  R  r  =  X  '  .  M  =•  > 

1010  IF  P  r  >  =  1  E  -  1  2  THE  N  Y  =  LGT  •:  F'  r  > 

1020  IF  P r <  = - 1 E - 1 2  THEN  Y=-24-LGT < -Pr > 

1030  IF  flES(Pr)< IE- 12  THEN  Y=-12 
1040  PLOT  fli,  r’ 

1050  NEXT  Ms. 

1060  PENUP 

1070  DUMP  GRAPHICS 

1080  OUTPUT  0;"  11 

1090  IF  ln>0  THEN  1200 

1100  FOR  M £  =  M 2  TO  Ml 

1110  IF  Y  <  fl  s  )  <  =  0  THEN  1130 

1120  NEXT  Ms 

1 130  M  3  =  M  s - 1 

1140  REDIM  X-M2:M3> 

1150  FOR  M •=.  =  M 2  TO  M3 

116  0  X  <  M  s  >  =  F  N  I  n  v  p  h  i  (  Y  (.  M  s.  >  > 

1170  NEXT  Ms 

1  1 8  0  PRINT  #  1 ;  X  <  *  >  !  S  t  o  r  e  f  al  se  a  1  ar  m  p  r  o  b  ab  i  1  i  t  y 

1190  GOTO  1380 
1200  REDIM  X<M2:M3> 

1210  READ  # i ; X ( * >  !  Read  i n  f a 1 se  al arm  probab i  1  i  t y 

1220  I d  = I d  e 1 1 a  1 1 n  > 

1230  J2=M2/Id 
1240  J 3  = I N T  < M 3 / I d  >  +  1 
1250  FOR  J = J 2  TO  J3 
1260  Y  <  J  )  =  F  M I  n  v  p  h  i  C  Y  <  J  )  > 

1270  NEXT  J 

1280  PLOTTER  IS  "9372A" 

1290  LIMIT  X 1 , X 2 , Y 1  ,  Y  2 
1300  OUTPUT  705; -'VSS'1 
1310  SCALE  S , O , S , E 
1320  FOR  J  =  J 2  TO  J3 
1330  T = J  *  I  d 
1340  IF  T > M 3  THEN  1370 
1350  PLOT  X  <  T  )  ,  Y  <  J  ) 

1360  NEXT  J 
1370  PENUP 
1380  NEXT  In 
1390  END 


1400 

1410 

1460 

SUE 

Exp X ,  Y ,  A  ,  E  ) 

!  EXPCZ) 

1470 

SUE 

L  o  q  <  X ,  Y ,  fl ,  B  ) 

!  PRINCIPAL 

1550 

1  560 

! 

DEF 

FNInvphi <X) 

!  INVPHI 

<  X  > 

via  AMS  55, 

1570 

IF 

X > 0  )  AND  <X<1) 

THEN  1600 

1580 

P  =  9 

. 99999999999E90* 

C  2  *  X  - 1 ) 

1590 

GOTO  1670 

1600 

IF 

X  = . 5  THEN  RETURN 

0 

1610 

P  =  X 

1620 

IF 

X > . 5  THEN  P=.5-< 

X- .  5« 

1630 

P  =  S 

Q  R  <  -  2  *  L  0  G  <  P  )  > 

1  640 

T=l+P*< 1 . 432788+P*< 

. 1 89269+P* . 

0  0  1 

3  0  8  >  > 

1650  P  =  P - <  2 . 515517+P*< . 802853+P*. 010328) >xT 

1660  IF  X< . 5  THEN  P  =  -P 

1670  RETURN  P 

1680  FNEND 

1690  ! 

1  7  0  0  S  U  B  F  f*  1 1  3  2  (  N  ,  X  *  *• )  ,  Y  K  *  )  )  1  N  <  =  2  *"•  1  ’J  , 
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LOG  <  Z ) 
26.2.23 


N  =  2  INTEGER,  O  S  U  B  S  C  R I  P  T 
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APPENDIX  E.  ASYMPTOTIC  EXPANSIONS  FOR  DISTRIBUTIONS  WHEN  r> 0 


The  characteristic  function  of  interest  is  given  by  (100)  and  (102): 


fh(f)  =  (l+if)  (l-1fu»)  exp 


T-^j  ’ 


(E-l) 


where  for  notational  convenience  in  this  appendix,  we  let 


v  =  ^j-,  n  =  Nr^,  o)  =  1+2R  .  (E-2) 

The  cumulative  distribution  function  is  obtained  by  substitution  of  (E-l)  in 
(63): 


Ph(u)  - 


-1 

727 


a-l 


d r  r  (i+->i) 


v-V 


exp 


Mil 


-  iuf] 


The  v-th  powers  are  principal  value,  being  positive  real  where  C+ 
the  positive  imaginary  axis. 


crosses 


( E— 3 ) 


Now  let  z  =  l+i£  in  (E-3),  yielding 

Ph(u)  =  777  J  dz  Z_V  d+—)-V  exp[feHl  -  U(Z-1H  •  (E-4) 

C1 

The  contours  C+  and  in  (E-3)  and  (E-4)  are  depicted  as  dashed  lines  in 
figure  E-l.  The  pole  at  J=0  is  moved  to  z=l;  the  remaining  singularities  are 
branch  points  (v  non-integer);  the  v-th  powers  are  positive  real  where 
crosses  the  positive  real  axis.  For  u<0,  an  equivalent  contour  to  C-^  is 
that  indicated  by  C£  in  figure  E-l,  since  the  exponential  in  (E-4)  furnishes 
rapid  decay  in  the  left-half  z-plane.  We  write  (E-4)  in  the  form 

Ph(u)  =  727^  f  dz  Z"V  exP("uz)  9i ( z )  ,  ( E— 5 ) 


where 


E-l 
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9l(z)  =  (1-z)"1  ( l+iu— tuz )— v  exp[^=^]  . 


Figure  E-l.  Contours  of  Integration  for  Cumulative  Distribution  Function 

In  order  to  get  the  asymptotic  development  of  (E-5),  we  expand  g^(z)  in 
a  power  series  in  z. 


g^z)  =  g^O)  +  g{ (o)  z  +  ...  , 


(E-7) 


where 


9i  (0) 


/  \  yi  \ u  / 

gx(0)  -  (l+»rv  exp(jjl),  =  1  +  ♦  -y-g 


Appeal  to  [10,  p.  96,  (4)]  then  yields  (for  all  v) 


( E— 8) 


Ph{u)  '  «P<u>|j^f-  gx(0)  +  g;(0)J  = 

-  Hv)-1  (l+u)-v  (-P)'"1  exp[u  -  [l  -  4(l  ♦  ft  * 


(1+0,)' 

as  u  -  °°  .  ( E— 9 ) 


Substitution  of  (E-2)  in  (E-9)  then  yields  result  (104) 


E-2 
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When  (E-l)  is  substituted  into  (64)  instead,  we  obtain  the  exceedance 
distribution  function  in  the  form 


1  "  Ph(u)  =  T27  J  ^  5  1  (1+i^~V  exP 


[rife  -  H 


(E-10) 


Now  let  z  =  -  if,  to  get 


1  -Ph( 


u)  -TJT  5 


dz  (z  -I)  1  (1  +^-z)_V  (»z)  V  * 


exp 


♦  u  ( z  -  i) 

2  a) 

L  gj  z  -J 


( E— 11 ) 


The  contours  C_  and  C3  in  (E-10)  and  ( E— 11 )  are  depicted  as  dashed  lines 
in  figure  E-2.  The  pole  at  f=0  is  moved  to  z  =  1  /to ;  the  remaining 
singularities  are  branch  points. 


Figure  E-2.  Contours  of  Integration  for  Exceedance  Distribution  Function 


For  u>0,  an  equivalent  contour  to  C3  is  that  indicated  by  C4  in 
figure  E-2,  since  the  exponential  in  ( E— 1 1 )  furnishes  rapid  decay  in  the 
left-half  z-plane.  We  write  ( E— 1 1 )  in  the  form 


E-3 
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1-Ph(u)  .  .  (1M-V  exp  TJT  {  dz 


o)  Z 


where 


g2U)  -  (l  - 


( E— 13) 


As  above,  we  expand  g£(z)  in  a  power  series  in  z. 


( E— 14 ) 


and  substitute  it  in  (E-12);  employment  of  [10,  p.  105,  (2)]  now  yields  (for 
all  v) 


( E— 15 ) 


as  u  -»  +  °o. 


where  Iu(z)  is  the  modified  Bessel  function  of  the  first  kind.  This  is  the 
general  result  for  the  exceedance  distribution  function;  the  various 
parameters  given  in  (E-2)  relate  it  back  to  the  problem  of  interest  in  the 
main  text. 

As  a  check  on  this  result,  we  let  r-*0;  then  n-*0,  and  ( E— 15 )  reduces, 
via  [6,  9.6.7],  to 


as  u->+<a<>;  r=0.  ( E— 16 ) 

Employing  the  identifications  in  (E-2),  there  follows  from  ( E— 16 ) 


E-4 


TR  7045 


N-2 


2  ( 1 +R )  p(  N  /  2 ) 


1+2R 


*  1  +  —  ...  Ill  tJ_\H 

r  1+R  8u 


1+2R  (N-2) (N+4+4R) 


as  u  ;  r=0  . 


(E-17) 


In  order  to  compare  this  result  with  that  for  y=1,  sample  mean  removal,  we 
must  replace  N  here  by  N-l;  see  the  paragraph  under  (101).  When  this  is  done, 
(E-17)  reverts  precisely  to  (61)  and  (62). 

Returning  to  the  general  result  for  the  exceedance  distribution  function 
in  (E-15),  if  we  keep  n>0  and  use  [6,  9.7.1]  for  large  arguments  of  Iy(z), 
there  follows  the  simpler  (less  accurate)  result 


( E— 18 ) 


as  u-» 


When  (E-2)  is  substituted  in  (E-18),  the  result  quoted  in  (105)  follows. 
As  a  special  case  of  (E-18),  if  v=l  (i.e.  N=2),  then 


[- 


as  u  — ►  +°*  ;  v=l.  ( E— 19 ) 


E-5/E-6 
Reverse  Blank 
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APPENDIX  F.  EXCEEDANCE  DISTRIBUTION  FUNCTION  FOR  Y=0,  N=l,  r>0 
CHARACTERISTIC  FUNCTION  APPROACH 


When  characteristic  function  (100)  with  N=1  is  substituted  in  (64),  the 
expression  for  the  exceedance  distribution  function  becomes 


1  - 


Ph(u) 


'  *2"  ■[ 


o-l 


-1/2 


d?  p  (1+i^)  d-iJo>) 


*-1/2 


exp 


1  fr 


l-ijo 


i 


-  iuf  ,  (F-l) 


where  w  =  1+2R  as  in  (102).  The  square  roots  are  taken  as  +1  at  |=0.  For 
u>0,  the  contour  C_  can  be  modified  to  that  indicated  in  figure  F-l,  where 
the  contributions  of  the  large  circular  arcs  in  the  lower-half  f-plane  tend  to 
zero.  The  small  circle  of  radius  p  centered  at  branch  point  -i/w  must  have 
0  <  p  <  l/u),  since  the  latter  is  the  distance  to  the  pole  at  the  origin. 


Figure  F-l.  Equivalent  Contours  for  (F-l) 


It  is  easy  to  show  that  the  two  vertical  contributions  in  figure  F-l  are 
equal.  Under  the  change  of  variable 


(F-2) 


F-l 
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the  sum  of  the  two  vertical  contributions  to  exceedance  distribution  function 
(F-l)  becomes 


+co 


(F-3) 


This  integral  remains  convergent  even  as  p-*0+.  Furthermore,  the  integrand 
decays  rapidly,  has  no  cusps,  and  involves  only  elementary  functions;  also  the 
integral  is  a  sum  of  positive  quantities  and  retains  significance  even  for 
large  u. 

On  the  small  circular  contour  Cp  in  figure  F-l,  let 


( F— 4 ) 


to  obtain,  for  the  circular  contribution  to  the  exceedance  distribution 
function,  the  quantity 


IT 


( F— 5 ) 


where  we  define  in  this  appendix 


(F-6) 


a) 


F-2 
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The  exceedance  distribution  function  is  given  exactly  by  the  sum  of  (F-3) 
and  (F-5),  for  any  0  <  p  <  1  /to .  It  would  be  advantageous  numerically  to  let 
p-»0+  in  these  two  equations;  however,  the  limit  of  (F-5)  is  not  obvious  and 
can  easily  be  done  incorrectly. 

AN  ERRONEOUS  APPROACH  FOR  C(p) 

It  is  tempting  to  let  p-*0+  in  those  locations  in  (F-5)  where  it  will  "do 
no  damage",  obtaining  for  the  integral  with  scale  factor  p^'^  the  quantity 


TT 


-7T 


(The  fallacy  of  doing  this  for  a  residue  calculation  with  an  essential 
singularity  is  demonstrated  in  the  next  subsection.)  Furthermore,  the  limit 
of  (F-7)  as  p-»0+  can  in  fact  be  determined  in  closed  form,  as  follows. 

Observe  that  the  integrand  of  (F-7)  has  a  saddle  point  in  the  complex  9-plane 
at 


& 


s 


=  -i  L  , 


(F-8) 


this  is  in  fact  the  only  saddle  point  in  the  (  —it, it )  strip  in  the  9-plane.  Now 
let  z  =  9-9$,  getting  for  (F-7) 


I 

p 


ir+iL 


-ir+iL 


dz  expjj  j  +  j  exp(-iz)J 


(F-9) 


The  radius  p  now  appears  only  in  the  limit  L  of  the  integral,  and  the 
integrand  has  a  saddle  point  at  z=0. 


F-3 
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The  straight  line  contour  for  (F-9)  can  be  deformed  into  contour  C, 
depicted  in  figure  F-2,  which  goes  through  the  saddle  point  at  z=0.  Now  if 


o-*-0+,  then  +  and  (F-9)  yields 


TT  +  ioD 


In  =  (20 


1/2 


dz  exp 


[i  f  +  \  exp(-iz)J  ,  ( F— 10) 


-TT+IOO 


where  the  contour  is  the  limit  of  C  in  figure  F-2  as  that  is,  the 

integral  is  between  the  two  valleys  at  ±ir+ico  and  is  connected  by  the  saddle 
point  at  z=0. 

The  steepest  descent  curves  out  of  the  saddle  point  of  the  integrand  of 
(F-10)  are  given  explicitly  by 

y  =  _in(-^-^)  for  -it  <  x  <  *  .  ( F— 1 1 ) 


Thus  if  we  let 


F-4 


TR  7045 


z=x  +  iy  =  x-  i  ln(s1  v-x-) 


&l_  _  ,  .  /Cos  x  1» 

dx  -  'sin  x  "  x'  ’ 


(F-12) 


on  the  steepest  descent  curves  in  (F-10),  there  follows 


I0  =  (2x)1/2 


COS  X  + 

sin  x 


+ 


1  x 
7  sin  x 


exp(-ix 


3 


-7T 


.  COS  X  + 

sin  x 


TT 


-  (2x)1/2  2  jdx  (Sii-i)1'2  exp 
0 

(The  integral  value  of  (2ir)^2 


(r sin  x)  "  (2x)1/2  2(2,r)1/?  =  4(\ir)1  /2  .  ( F— 13 ) 

in  ( F— 13 )  was  deduced  by  numerical  integration.) 


Recalling  the  definition  of  I  in  (F-7),  we  then  have  the  dubious 
result  for  the  limit  of  (F-5): 


(1+R) 


1/2 


( F— 14) 


where  we  employed  (102)  and  (F-6).  Actual  numerical  evaluation  of  (F-14), 
combined  with  V(0)  from  (F-3),  gives  incorrect  results  for  the  exceedance 
distribution  function  ( F— 1 ) ;  thus  the  replacement  of  p  with  0  in  (F-5)  is 
invalid.  The  explanation  for  this  pitfall  is  the  essential  singularity  of 
(F-l)  at  £  =  — i / o> ;  a  simpler  illustration  follows. 


F-5 
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RESIDUE  OF  ESSENTIAL  SINGULARITY 


The  function 


exp(j) 


( F— 15 ) 


has  an  essential  singularity  at  z=0,  with  residue  1,  as  exemplified  by  this 
Laurent  expansion.  Now  consider  the  function 


f ( z )  =  exp(i)  g(z)  , 


(F-16) 


where  g(z)  is  analytic  at  z=0.  Then 


f (z )  =  (i  +  7  +  ^77  +  •••  Yg(o)  +  g(1)(o)  z  +  t^(2)(o)  z2  +  ...j  . 
The  coefficient  of  1/z  in  ( F— 17 )  is  the  residue  of  f(z)  at  z=0;  namely 


(F-17) 


Res  =  g ( 0 )  +  jr  9^(0) 


+O0 


-5 

n=0 


9(n)(0) 

ni  (n+1) 1 


(F-18) 


Thus  the  residue  of  f(z)  at  z=0  depends  on  the  behavior  of  g(z)  in  a 
neighborhood  of  z=0,  and  not  just  the  value  g(0). 

A  couple  of  examples  yield  the  following: 

g(z)  =  (l-az)-1  ,  Res.SiEiiH.; 

2  Ii(2a) 

g(z)  =  exp(a  z)  ,  Res  = - - —  . 


(F-19) 


F-6 
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CORRECT  APPROACH  FOR  C(p) 

Reconsider  the  integral  in  (F-5)  plus  the  scale  factor  p^2;  making  the 
substitution  z  =  ©+iL,  where  L  is  given  in  (F-8),  there  follows  for  this  quantity 


n+iL 


-TT+i  L 


dz  (l-2wx  eiz)-1  (1 


2m\ 

1+0) 


*  exp[i 


1 

2 


+  2xu  eiz]  . 


( F— 20 ) 


The  uppermost  singularity  of  the  integrand  in  the  z-plane  (within  the  -it,  ir 
strip)  is  a  pole  at  Zp  =  i  Jn(2wx);  however,  the  straight  line  contour  in 
(F-20)  remains  above  this  pole  because  p  <  1/w;  see  (F-8).  Furthermore,  the 
total  integrand  of  (F-20)  has  a  saddle  point  on  the  imaginary  axis  of  the 
z-plane  above  the  pole  location  zp,  because  the  integrand  is  infinite  at  the 
pole  and  at  z  =  0  +i».  Thus  the  straight  line  contour  in  (F-20)  can  be 
modified  so  as  to  pass  through  the  saddle  point,  and  yet  remain  above  Zp. 
Fina^y,  letting  p-»0+,  then  L-^+o*,  and  (F-20)  combined  with  (F-5)  yields  the 
exact  result  for  the  circular  component 


C(0) 


n+ioo 


-ir+i<» 


(l-2u)X  e1z)  * 


*(X  “  vt  e1j  U2  exp£  f  +  \  e_iZ  +  2xu  eiZJ  •  (p-21) 

where  the  two  valleys  of  the  integrand  at  ±ir+i«o  are  joined  with  a  contour 
through  the  saddle  point  lying  above  the  pole  at  zp  _  j  j^n(2aix).  Here 
a)  =  1+2R. 


F-7 
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The  other  component  of  the  exceedance  distribution  function, 
corresponding  to  (F-21),  is  given  by  (F-3)  at  p  =  0+: 


0 


dt  (1+t2)  * 


(F— 22) 


Thus  for  u  >_  0,  (F-l)  and  figure  F-l  yield  exceedance  distribution  function 


1  -  Ph(u)  =  C(0)  +  V(0)  =  (F-21)  +  (F-22)  . 


( F— 23 ) 


Computationally,  (F-21)  is  not  too  attractive,  because  of  the  complex 
integrand  and/or  the  need  to  determine  the  steepest  descent  paths  to  ^ir+i0* 
numerically.  Accordingly,  an  alternative  direct  procedure  for  determining  the 
exceedance  distribution  function  of  random  variable  h  is  now  presented. 

DIRECT  EVALUATION  OF  EXCEEDANCE  DISTRIBUTION  FUNCTION 

For  y=0,  N=l,  (12)  and  (3)  yield  the  crosscorrel ator  output  for  the 
signal  and  noise  model  as 


q  =  u1v1  =  [uu  +  us(l)  +  ud(l)]  [uv  +  vs(l)  +  vd(l)]  ..  ( F— 24 ) 


The  normalized  crosscorrelator  output  is  then,  from  (49)  and  (97), 


where  x  and  y  are  joint  Gaussian  random  variables  witn  statistics 


F-8 
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x  =  r  ,  y  =  r  , 
u  ’  J  v  ’ 


°X  -  1  +  Ru’  «y  -  1  +  V  (x-x)(y-y)  -  PS(RURV)1/Z  • 


We  now  make  the  same  assumptions  as  in  (99);  see  also  (56)  et  seq. 
(F-26)  specializes  to 


x  -  y  -  r,  a2  “  ay  ”  1  +  R  r  °2’  pxy  “  1+R  ’ 

The  joint  probability  density  function  of  x,  y  is  then  given  by 


P2(x, 


Ul-.  2)1/fl  1  exp 

(x-r)2  +  (y-r)2  -  2pxy(x-r) (y-rfj 

ld"  l,1  °xy/  J  exp 

2«2(1  -  px2)  J 

We  now  have  cumulative  distribution  function 


Ph(u)  =  ^  dx  dy  p2(x,y)  for  u  <  0  , 


r2+r4 


and  exceedance  distribution  function 


1  -  Ph(u)  =  f j'  dx  dy  p2(x,y)  for  u  >  0  , 


^l+^3 


where  regions  Rj,  r2,  r3,  r4  are  indicated  in  figure  F-3. 


(F-26) 

Then 


(F-27) 


(F-28) 


(F-29) 


(F-30) 


F-9 
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Figure  F-3.  Regions  of  Integration 


If  we  now  rotate  axes  according  to 


s 


(P-31) 


and  employ  (84)  and  (F-27),  there  follows  (after  scale  changes  of  the 
vari ables) 


Ph(u)  =  2  |  dv  0(v  -  r(2/u) i;^)  $^(u>v2-2u)1/2J 


for  u  <  0 


( F— 32 ) 


(F-33) 


Here  u  =  1+2R.  These  real  integrals  are  very  useful  for  the  evaluation  of  the 
distributions  of  h  when  y=0>  N*l.  In  fact,  (F-33)  is  preferred  over 
(F— 21 )— (F— 23) ;  but  (106)  is  preferred  over  (F-32)  since  $  need  not  be 
evaluated  in  (106).  This  is  in  fact  the  procedure  utilized  here  to  obtain 
numerical  results  for  this  case  of  y=0»  N=l. 
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APPENDIX  G.  PROGRAM  FOR  EVALUATION  OF  OPERATING  CHARACTERISTICS  FOR  y=0 

The  comments  in  appendices  C  and  D  are  relevant  here  also.  The 
characteristic  function  used  as  the  starting  point  is  given  b.y  (100).  It  was 
observed  under  (103)  that  jT^(f)|  for  (100)  is  monotonical ly  decreasing  for 
all  $>  0;  thus  the  choice  of  truncation  value  L  is  simplified;  see  appendix  D 
comments.  A  table  for  sampling  increment  aq  (when  R=0)  follows. 


N 

aq  for  r=l 

> 

O 

O 

'S 

3 

.07 

o 

o 

(_n 

4 

.07 

.05 

8 

.05 

.03 

16 

.04 

.02 

32 

.03 

.02 

64 

.025 

.015 

128 

.020 

.010 

256 

.012 

.007 

Table  G-l.  Values  of  aq  for  y=0 


10  ' 

GAMMA  =  0  NO 

SAMPLE  MEAN  REMOVAL 

20 

He  =32 

| 

N ,  N  u  fn  b  e  r  o  f  t  e  r  rn  s  a  d  d  e  d 

30 

R  s  =  1 

i 

r ,  N o  r  rn  a  1  i  zed  rn e  an 

4  0 

D  e  1 t  a0  =  .  03 

i 

I  r'i  i  t  i  a  1  d  e  1 1  a 

50 

E •=•  =  2 * P  I  /He  1  t.  a0 *  .  37 

5  ! 

Bias  b  <  d  e  p  e  n  d  s-  o  n  r  ) 

60 

M  f  =  2  A  1 0 

j 

Size  of  FFT 

70 

OUTPUT  0 ; "GAMMA  = 

0 " ; "  N 

=  " ; He ;  "  r  = " ; Rs 

80 

OUTPUT  0 ; " “ 

90 

DATA  -4, -3, -2. 5,-2 

,  — 1-5,— 1,— 

.5,0, .5,1, 1.5,2 

100 

READ  Ns  <  *  > 

i 

SNR  R=2An 

110 

output  0 ; Us <  *  > ; 

120 

DATA  1,2, 2, 2,2, 2, 2 

,  2 , 2 , 4 , 4 , 4 

130 

READ  Idelta<*> 

140 

M  AT  Del  t  a=  < D e  1  t  a0  > 

/•  I  de  1 1  a 

150 

OUTPUT  0;  BeHaC*) ; 

160 

DATA  IE- 10, IE-9,  IE 

-8, 1 E-7 , IE 

- 6 , 1 E — 5 , 1 E — 4 , .001, .01,  . 1 ,  . 5 , . 9 , 

170 

READ  S c (* ) 

180 

DIM  Hs< 1 : 12) , I del t 

a<0: 12) , De 

1 t  a<  0 : 12), Sc  Cl: 1 4 > 

19U 

D I M  X <0:319l>, V < 0 : 

8  191) 

200 

FOR  1=1  TO  14 

210 

S c < I )  =  F HI n m p h i <  S c  < 

I  )  ) 

220 

NEXT  I 
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230 

S=Sc  l  > 

2  4  0 

E  =  S  c  (  1  4  > 

250 

S'  ■:  a  1  e  =  C  E  -  S  >  x  <  0  -  S  ) 

260 

X  1=30 

270 

X  2  =  1  7  0 

280 

VI  =35 

290 

Y2=Y  1  +  (  X2-X 1  )  * S ■:  al  e 

300 

PLOTTER  IS  " 9872 H “ 

3 1  0 

LIMIT  X 1 5 X  2 , Y 1 , Y  2 

320 

OUTPUT  705; “VS 3" 

330 

SCALE  S , 0  , S , E 

340 

FOR  1=1  I 0  14 

350 

MOVE  S,Sc<n 

360 

DRAW  0 , Sc  < I > 

370 

NEXT  I 

380 

FOR  1=1  TO  11 

390 

MOVE  Sc<I),S 

400 

DRAM  5c  (.  I  )  ,  E 

410 

NEXT  I 

420 

MOVE  $,S 

430 

DRAW  0,8 

440 

PEHUP 

450 

M  1  =  M  f  - 1 

466 

N2=Nc/2 

470 

R  s  n  =  N  c  *  R  s  *  R  s 

N  r*  2 

480 

FOR  I n  =  0  TO  12 

490 

IF  I n > 0  THEN  520 

500 

Rc  =0 

510 

GOTO  530 

520 

Rc  *2ANs < I n )  1 

SNR  R  =  2 n 

530 

OUTPUT  0;“R  =“ ; Rc , “  Celt 

a  = “ ; D  cl t a< I n ) 

540 

ASSIGN  #1  TO  “ ABSC I S “  1 

T  e  rri  p  o  r  a r  y  s  t  o  r  ag  e 

550 

Belt  a= D e 1 t  a < I n >  ! 

f  o  r  f  a  1  «•  e  a  1  ar  rn 

probab i 1 i t y 

560 

R  2  =  R  c  *  2 

570 

R21=R2+1 

530 

Mux=Nc *Rc  +  Rsn  1 

M  e  a  n  o  t  r  a  n  id  o  m 

variable  h 

590 

Muy=Mux+Bs  ! 

M  e  an  o  f  s-  h  i  f  t  e  d 

variable,  y 

6  0  0 

RED  I M  X<0: Ml ) , Y<0S Ml > 

610 

MAT  X  =  Z E R 

620 

MAT  Y=ZER 

630 

X  < 0 ) =0 

6  4  0 

Y  <  0  >  = .  5  *  Ii  e  1 1  a*  M  u  y 

650 

Ls=0 

660 

Lft-Ls+1 

670 

X  i  =  Be  1 1,  a*Ls  ! 

A  r  g  u  m  e  n  t  x  i  o  f 

char.  fn. 

680 

E i =X i *R2 1 

j 

Cal  c u 1  at i on 

690 

CALL  Log < 1 +X i ^E i , -X i *R2 , A i , 

Ei  )  ! 

o  f 

700 

CALL  D i v < 0 , X i *Rsn , 1 , -E i , C i , 

Di  >  ! 

c  h  ar  ac  t  e  r  l  s 

710 

CALL  Exp*:  Ci  - N 2  *  A  i  ,  Di  +Xi  *Es- 

N2*B i , Fyr , Fy i >  ! 

f  u  n  c  t  i  o  n 

720 

M  £•  =  L  £■  MOD  M  f 

{ 

fy<xi ) 

730 

Ar=FyrxLs 

740 

Ai =Fy  l xLs 

750 

X  C  M  s  >  =  X  <  M  s  )  +  A  r 

760 

Y  c;  M  s  >  =  Y  (.  M  s  )  +  A  i 

770 

M  ag  s  q  =  A r  *  H  r  +  A  i  * H  i 

780 

I F  *  M ag s q  > 1 E - 2 4  THEN  660 

790 

OUTPUT  0 ;  “  X  i  = 11 ;  X  i  ;  11  Mag 

=  "  ;  SQR  *•  Magsq  > 

800 

CALL  F  ft  1 3z <Mf , X  <  *  > , Y  C * >  > 
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810  FOR  M s - 0  TO  Ml 
8 20  T  =  V < M  s >  • '  P I - M s ✓ M f 

8  3  0  X  (  M  s  >  = .  5  -  T  !  C  u  rn  u  1  at  i  u  e  d  i  s  t  r  i  b  u  t  i  o  n  f  u  n  c  t  i  o  n 

8  4  0  V  M  s  >  =  .  5  +  T  !  E  x  c  e  e  d  an  c  e  d  i  s  t.  r  i  b  u  t  i  o  n  f  u  n  c  t  i  o  n 

850  NEXT  Ms 

860  OUTPUT  0 ; Y < 0 > ; Y <  1 ) ; Y < M 1 - 1 > ; Y < M 1 ) 

870  PLOTTER  IS  "GRAPHICS” 

880  GRAPHICS 
890  SCALE  0«Mf,-!4,0 
900  LINE  TYPE  3 
910  GRID  M f / 8 ,  1 
920  PENUP 
930  LINE  TYPE  1 
940  FOR  M s = 0  TO  Ml 
950  Pr=Y < Ms  > 

960  IF  P  r  >  = 1 E - 1 2  THEN  Y  =  LGT  <Pr  > 

970  IF  P  r  <  =  -  1  E  -  1  2  THEN  Y=-24-LGT < -Pr ) 

980  IF  ABSXPrX  IE- 12  THEN  Y  =  -12 

990  PLOT  Ms , Y 

1000  NEXT  Ms 

1010  PENUP 

1020  FOR  Ms  =  0  TO  Ml 

1030  P  r  =  X < M s ) 

1040  IF  P r >  = 1 E - 1 2  THEN  Y  =  LGT <Pr ) 

1050  IF  P r <  =  -  1  E -  1  2  THEN  Y— 24-LGT<-Pr ) 

1060  IF  ABS<* Pr^lE-12  THEN  Y  =  -12 

1070  PLOT  M s , Y 

1080  NEXT  Ms 

1090  PENUP 

1160  DUMP  GRAPHICS 

1110  OUTPUT  0;  "  " 

1120  IF  I n > 0  THEN  1270 

1130  FOR  M s  =  0  TO  Ml 

1140  IF  Y ( M s  > < . 7  THEN  1160 

1150  NEXT  Ms- 

1160  M  2  =  M  s  - 1 

1170  FOR  M s  =  M 2  TO  Ml 

1180  IF  Y < M s > < =  8  THEN  1200 

1190  NEXT  Ms 

1 200  M 3  =  M  s - 1 

1210  REDIM  XCM2:M3> 

1220  FOR  Ms=M2  TO  M3 
1230  X < M  s > =  F N I n  m p h i < Y  < M  s  > > 

1240  NEXT  Ms 

1250  PRINT  #  1  ;  X  (  *  >  !  S  t  o  r  •=  f  a  1  s  ■=•  a  1  a r  rn  p  r  obab  i  1  i  t  y 

1260  GOTO  1460 
1270  REDIM  X<M2:M3> 

1280  READ#1;X<*>  1  Read  in  false  al a r m  p r o b ab i 1 i t y 

1290  I  d  =  I  d  e  1 1  a  <  I  n ') 

1300  J 2  = I N T  <  M 2 x I d  > 

1310  J  3  = I N  T ( M  3  x I d  >  +  1 
1320  FOR  J=J2  TO  J3 
1330  Y  <  J  >  =  F N I n m p h i <  Y  < J > > 

1340  NEXT  J 
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1350 

PLOTTER  IS  "  93721=1 " 

136© 

LIMIT  X1,X2,Y1,Y2 

1370 

OUTPUT  705; " VS3" 

1380 

SCALE  S , 0 , S ,  E 

1390 

FOR  J  =  J2  TO  J 3 

1400 

T  -  J  *  I  d 

1410 

IF  T < M 2  THEM  1440 

1420 

IF  T > M 3  THEM  1450 

1430 

PLOT  X < T > , V < J > 

1440 

NEXT  J 

1450 

PENUP 

1460 

NEXT  In 

1470 

END 

1480 

i 

1490 

SUE  D i v < X 1 , V 1 , X2 , Y2 , A , B > 

!  D I  V  <  Z  > 

1540 

! 

1550 

SUE  Exp  < X , Y ,  A ,  E  > 

•  E  X  P  <  Z  > 

1600 

i 

1610 

SUE  L o g ( X , Y, A, B> 

'  PRINCIPAL  LOG<Z> 

1690 

! 

1700 

D  E  F  F  N  I  n  v  p  h  i  <  X  > 

'  INVPHI <X>  via  AMS  55,  26.2.23 

1830 

i 

1840 

SUE  Ff 1 1 3 z  < N , X < * ) , Y < *  >  > 

!  N  <=  2713,  N=2A INTEGER,  0  SUBSCRIPT 

f 

V 

G-4 


TR  7045 


REFERENCES 


1. 

A.  H.  Nuttall,  "Resolution  of  Ambiguity  for  Randomly  Moving  Line  Array," 

NUSC  Technical  Memorandum  831150,  3  October  1983. 

2. 

T 

A.  H.  Nuttall,  "Exact  Performance  of  General  Second-Order  Processors  for 
Gaussian  Inputs,"  NUSC  Technical  Report  7035,  15  October  1983. 

3. 

A.  H.  Nuttall,  "Accurate  Efficient  Evaluation  of  Cumulative  or  Exceedance 
Probability  Distributions  Directly  from  Characteristic  Functions,"  NUSC 
Technical  Report  7023,  1  October  1983. 

4. 

A.  H.  Nuttall  and  B.  Dedreux,  "Exact  Operating  Characteristics  for  Linear 
Sum  of  Envelopes  of  Narrowband  Gaussian  Process  and  Sinewave,"  NUSC 
Technical  Report  7117,  11  January  1984. 

5. 

I.  S.  Gradshteyn  and  I.  M.  Ryzhik,  Table  of  Integrals,  Series,  and 
Products,  Academic  Press,  Inc.,  New  York,  1980. 

6. 

Handbook  of  Mathematical  Functions,  U.  S.  Department  of  Commerce, 

National  Bureau  of  Standards,  Applied  Mathematics  Series  No.  55, 

U.  S.  Government  Printing  Office,  June  1964. 

7. 

A.  H.  Nuttall,  "Some  Integrals  Involving  the  Q-Function,"  NUSC  Technical 
Report  4297,  17  April  1972. 

'  8. 

A.  H.  Nuttall,  "Some  Integrals  Involving  the  QM-Func.tion, "  NUSC 

Technical  Report  4755,  15  May  1974. 

9. 

W.  B.  Davenport  and  W.  L.  Root,  Random  Signals  and  Noise,  McGraw-Hill 

Book  Co.,  N.Y.,  1958. 

10. 

N.  W.  McLachlan,  Complex  Variable  Theory  and  Transform  Calculus, 

Cambridge  University  Press,  second  edition,  N.Y.,  1963. 

R-l/R-2 
Reverse  Blank 


TR  7045 


INITIAL  DISTRIBUTION  LIST 

Addressee  No.  of  Copies 

ASN  ( R&S)  1 

OUSDR&E  (Research  &  Advanced  Technology)  1 

Deputy  USDR&E  (Res  &  Adv  Tech)  1 

Principal  Dep  Assist  Secretary  (Research)  1 

OASN ,  Dep  Assist  Secretary  (Res  &  Applied  Space  Tech)  1 

ONR,  0NR-100,  -102,  -200,  -400,  -410,  -425AC,  -430  7 

CNO,  OP-098 ,  -951  ,  -951D,  -951E,  -96,.  -981F  6 

CNM,  MAT-00,  -05,  -05B  (CAPT  Z.  L.  Newcomb),  SP-20  4 

DIA,  DT-2C  1 

NRL  1 

NRL,  USRD  1 

NRL,  AESD  1 

NORDA  1 

•  USOC,  Code  241,  Code  240  (Acoustics  Research  8ranch)  2 

NWOC  1 

NAVOCEANO  1 

NAVOCEANO,  Code  02  (Scientific  and  Technical  Director)  1 

NAVELECSYSCOM,  ELEX  03  1 

NAVSEASYSCOM,  SEA-00  1 

NAVAIRDEVCEN,  Warminster  1 

NAVAIRDEVCEN,  Key  West  1 

NOSC,  Code  8302  1 

NOSC,  Code  6565  (Library)  1 

NAVWPNSCEN  1 

NCSC ,  Code  724  1 

NAVCIVENGRLAB  1 

NAVSWC  1 

NAVSURFWPNCEN,  Code  U31  1 

DWTNSRDC  ANNA  1 

DWTNSRDC  BETH,  Code  U31  1 

NISC  1 

CNET,  Code  017  1 

NAVTRAEQUIPCENT,  Technical  Library  1 

NAVPGSCOL  1 

NAVWARCOL  1 

APL/UW  SEATTLE  1 

ARL/PENN  STATE,  STATE  COLLEGE  1 

CENTER  FOR  NAVAL  ANALYSES  (ACQUISITION  UNIT)  1 

DTIC  12 

DARPA  1 

NOAA/ERL  1 

NATIONAL  RESEARCH  COUNCIL  1 

WEAPON  SYSTEM  EVALUATION  GROUP  1 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION  1 

ENGINEERING  SOCIETIES  LIB,  UNITED  ENGRG  CTR  1 

NATIONAL  INSTITUTE  OF  HEALTH  1 

ARL,  UNIV  OF  TEXAS  1 

MARINE  PHYSICAL  LAB,  SCRIPPS  1 

UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO  1 

NAVSURWEACTR  1 


TR  7045 


INITIAL  DISTRIBUTION  LIST  (Cont'd) 

Addressee  No.  of  Copies 

INTERMS  INC.  1 
MAR  INC.  1 
BBN  1 
EWASCTRI  1 
HYDROINC  (D.  Clark)  1 
ANALTECHNS  1 
TRA  CORP.  (0.  Wilkinson)  1 
BB&N  INC.  1 
GESY  (D.  Bates)  1 
ORI  CO.,  Inc.  (G.  Assard)  1 
ROYAL  MILITARY  COLLEGE,  Kingston,  Ontario,  Canada  (Prof.  Chan)  1 
UNIVERSITY  OF  FLORIDA,  Dept,  of  Elec.  Eng.,  (D.  Childers)  1 
ROCKWELL  INTERNATIONAL,  Anaheim,  CA  (Dr.  D.  Elliott)  1 
NRL,  Systems  Engineering  and  Advanced  Concepts  Branch 

(A.  Gerlach)  1 
BB&N,  Cambridge,  MA  (H.  Gish)  1 
Bernard  Harris,  Dobbs  Ferry,  NY  1 
NOSC ,  SAN  DIEGO,  Code  713  (F.  Harris)  1 
EG&G,  MANASSAS,  VA  (Dr.  J.  Hughen)  1 
DREXEL  UNIV.,  Dept,  of  Elec.  &  Comp.  Eng.,  Philadelphia, 

PA,  (Prof.  S.  Kesler)  1 
SCHLUMBERGER  WELL  SERVICES  ENGINEERING,  Houston,  TX 

(S.  Marple)  1 
Dr.  David  Middleton,  New  York,  NY  1 
BB&N,  Canoga  Park,  CA  (Dr  A.  Piersol)  1 
UNIV.  OF  RHODE  ISLAND,  Dept,  of  Elec.  Eng.  (Prof.  L.  Scharf, 

Prof.  D.  Tufts)  1 
YALE  UNIVERSITY,  Dept,  of  Elec.  Eng.,  New  Haven,  CT 

(Prof.  P.  Schultheiss)  1 
ARL,  PENN  STATE  UNIV.  (F.  Symons)  1 
RCA  Missile  and  Surface  Radar  (H.  Urkowitz)  1 
TELAVIV  UNIVERSITY,  Dept,  of  Electronic  Systems  (Prof. 

E.  Weinstein)  1 
NRL,  Washington,  DC,  Code  5162  (N.  Yen)  1 
UNIVERSITY  OF  IOWA,  Elec.  &  Comp.  Eng.  Dept.  (Prof.  D.  H.  Youn)  1 


